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Heat and Mass Transfer in Aeronautical 
Engineering 


D. B. SPALDING 
(Imperial College of Science and Technology) 


SUMMARY: The main purpose of the paper is to show that standard techniques 
of heat and mass transfer theory, as practised in other branches of engineering, 
can be helpfully applied to the corresponding aeronautical problems, particularly 
those arising in connection with high-speed flight and rocket propulsion. 


After a brief outline of the main features of the theory of heat and mass 
transfer in chemically reacting fluids, the standard techniques are applied to the 
following problems : — 


(a) The laminar boundary layer at the nose of an axi-symmetrical body at 
hypersonic flight speed. 


(() How to calculate the heat transfer rate to an internally cooled 
surface, taking account of dissociation. 


(ii) How to calculate the quantity of coolant needed for transpiration- 
cooling of the surface, taking account of possible combustion of the 
coolant. 


(ii) How to calculate the rate of burning of a graphite “heat-shield.” 
(iv) How to calculate the rate of “ablation” of a heat-shield material 
which melts and sublimes. 
(hb) The turbulent boundary layer in the nozzle of a rocket motor . 


(‘) How to calculate the heat flow from the gas to a bare metal nozzle 
wall:—simple formula; pipe-flow formula; boundary layer analysis. 

(i) How to calculate the rate of supply of liquid fuel to the motor wall 
if film-cooling is to be employed. 


(iii) How to estimate the thickness of carbon film which will be deposited 
on the wall of a liquid-oxygen-kerosine motor, and its influence on 
the heat transfer rate. 


(‘v) How to calculate the rate of enlargement of the throat area of a 
solid-propellant rocket motor with a graphite nozzle 


The procedures used comprise the following elements : — 
(a) Relations of the “Ohm’s Law” type, namely: 


“Mass transfer rate per unit area 
=surface conductance X dimensionless driving force 
=g xX B” 

and, in the absence of mass transfer, 


“Heat transfer rate per unit area=gX enthalpy dilference.” 


*A Lecture delivered to the Royal Aeronautical Society on 24th November 1959. 


(The Aeronautical Quarterly, Vol. XI, May 1960) 
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Heat and Mass Transfer in Aeronautical 
Engineering 


D. B. SPALDING 
(Imperial College of Science and Technology) 


SumMMaRY: The main purpose of the paper is to show that standard techniques 
of heat and mass transfer theory, as practised in other branches of engineering, 
can be helpfully applied to the corresponding aeronautical problems, particularly 
those arising in connection with high-speed flight and rocket propulsion. 


After a brief outline of the main features of the theory of heat and mass 
transfer in chemically reacting fluids, the standard techniques are applied to the 
following problems : — 


(a) The laminar boundary layer at the nose of an axi-symmetrical body at 
hypersonic flight speed. 
(i) How to calculate the heat transfer rate to an internally cooled 
surface, taking account of dissociation. 


(ii) How to calculate the quantity of coolant needed for transpiration- 
cooling of the surface, taking account of possible combustion of the 
coolant. 

(iii) How to calculate the rate of burning of a graphite “heat-shield.” 


(iv) How to calculate the rate of “ablation” of a heat-shield material 
which melts and sublimes. 


(b) The turbulent boundary layer in the nozzle of a rocket motor 


(i) How to calculate the heat flow from the gas to a bare metal nozzle 
wall:—simple formula; pipe-flow formula; boundary layer analysis. 

(ii) How to calculate the rate of supply of liquid fuel to the motor wall 
if film-cooling is to be employed. 

(iii) How to estimate the thickness of carbon film which will be deposited 
on the wall of a liquid-oxygen-kerosine motor, and its influence on 
the heat transfer rate. 

(iv) How to calculate the rate of enlargement of the throat area of a 
solid-propellant rocket motor with a graphite nozzle. 


The procedures used comprise the following elements : — 
(a) Relations of the “Ohm’s Law” type, namely : 


“Mass transfer rate per unit area 
=surface conductance X dimensionless driving force 
=gX B” 

and, in the absence of mass transfer, 


“Heat transfer rate per unit area=g enthalpy difference.” 


*A Lecture delivered to the Royal Aeronautical Society on 24th November 1959. 
{The Aeronautical Quarterly, Vol. XI, May 1960) 
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(b) Aerodynamic analysis leading to the evaluation of the conductance 
g, in Ib.mass/ft.*h. Boundary layer theory is used. 


(c) Thermodynamic analysis leading to the evaluation of B or of the 
enthalpy difference. Here diagrams of the enthalpy of the occurrent 
mixture versus its composition are used, permitting graphical evaluation 
of the required quantities. 


1. Introduction 
1.1. PURPOSE OF PAPER 


The terms “heat transfer” and “mass transfer,” as well as being the names of 
physical processes, also denote a body of engineering knowledge which is almost as 
extensive as the related subjects: thermodynamics and fluid dynamics. This 
knowledge has been built up during the last half-century, mainly by mechanical 
engineers concerned with boilers, heat exchangers, coolers and condensers for power 
plant, and by chemical engineers working on distillation, absorption, catalysis and 
other processes occurring in the chemical industry. 


In recent years the thermal problems of high-speed flight, and the practical 
measures used to overcome them, have brought aeronautical engineers to grips with 
problems of heat and mass transfer. Important advances have been made in this 
way. However, the available techniques are still only imperfectly exploited in 
aeronautical engineering, partly because of the differences of language and notation 
which divide one branch of engineering from another. 


The purpose of the present paper is two-fold: firstly, to demonstrate the 
existence of general and simple methods of calculating heat and mass transfer rates, 
which can deal with the aeronautical problems as special cases; and secondly, to 
show that these methods still need considerable development, which, in some 
respects, aerodynamicists are best fitted to provide. 


1.2. OUTLINE OF PAPER 


The points just mentioned will be illustrated by a discussion of two systems of 
current aeronautical interest. The first is the nose-cone of a high-speed missile; 
here the problem is to predict the rate of heat transfer from the dissociated gases 
adjacent the surface, and to assess the influences of transpiration cooling and of the 
so-called “ablation” of a “heat-shield.” 


The boundary layer on the nose-cone is laminar; the second system, by contrast, 
involves a turbulent boundary layer. The problem is the prediction of the heat 
transfer from the exhaust gases to the nozzle of a rocket motor. The mass transfer 
(erosion of the throat) which takes place when an uncooled graphite nozzle is used 
will also be discussed. 


It will, however, be necessary to preface discussion of these special problems by 
a short summary of some of the major relevant results of mass transfer theory. It is 
not possible to derive these in the present paper; reference will have to be made to 
the literature for justification of the general statements which will be used and for a 
discussion of the limits of their validity. 
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The final section lists some of the gaps in the current armoury of techniques 
for calculating mass transfer rates. 


NOTATION 
A 


b= 


B 


cross-sectional area of rocket nozzle (ft.?) 


(P—P;)/(Ps—Pr), dependent variable of the standard differential 
equation of mass transfer, equation (8) 


driving force for mass transfer appearing in “Ohm’s Law” of mass 
transfer, equation (1) 


drag coefficient 

specific heat at constant pressure (B.t.u./lb. mass °F) 
specific heat at constant volume (B.t.u./lb. mass °F) 
characteristic velocity of rocket propellant (ft. /sec.) 
diameter of rocket nozzle (ft.) 

function of B and (Pr) appearing in equation (21) 
Euler number 


function of B and (Pr) appearing in equation (21) 


mass fraction of a specified component in a mixture of two streams, 
irrespective of state of chemical aggregation 


surface conductance (lb. mass/ft.? h.) 
value of g for vanishing mass transfer (Ib. mass/ft.? h.) 


constant in Newton’s Second Law of Motion, 
=4-17 x 10° (Ib. mass ft./Ib. force h.”) 


mass velocity of main stream outside boundary iayer (Ib. mass/ft.’ h.) 
vectorial local total mass flux (Ib. mass/ft.? h.) 

point defining the stagnation condition of the combustion gases 
specific enthalpy of local mixture (B.t.u./1b. mass) (except in Fig. 4) 
local stagnation enthalpy (B.t.u./Ib. mass) 


enthalpy decrease of propellant gases where expanding isentropically to 
zero pressure (B.t.u./Ib. mass) 


mechanical equivalent of heat, = 778 (ft. Ib. force/B.t.u.) 

thermal conductivity of heat-shield material (B.t.u./ft. h. °F) 
mass fraction of chemical compound j in local mixture 

mass transfer rate into fluid under consideration (Ib. mass/ft.? h.) 


total rate of disappearance of heat-shield material through mass transfer 
and melting (Ib. mass/ ft.” h.) 


combustion chamber pressure (atm.) 


conserved property (e.g. A or f) 
107 
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Prandtl number of fluid stream 
Reynolds number 
distance of point on surface from axis of symmetry (ft.) 
radius of curvature of missile nose (ft.) 
heat flux by conduction to interior of missile nose (B.t.u./ft.? h.) 
heat flux across control surface L 
Stanton number, =2/G 
Stanton number for vanishing mass transfer, = g*/G 
absolute temperature (°F abs.) (except in Fig 4) 
temperature (°F) 
velocity of gas stream just outside boundary layer (ft. /h.) 

x distance along surface in the stream direction (ft.) 
distance normal to surface in direction of main stream (ft.) 


“exchange coefficient,” i.e. either diffusion coefficient multiplied by local 
mixture density or thermal conductivity divided by specific heat 
(Ib. mass /ft. h.) 


emissivity of radiating surface 

viscosity of gas (Ib. mass/ft. h.) 

density (Ib. mass /ft.*) 

Stefan’s constant, =0:1713 x 10-* [B.t.u./ft.? h. (°F abs.)*‘] 


main stream 

control surface on non-gas side of stream boundary 
coolant reservoir 

control surface on gas side of stream boundary 
transferred substance 


far upstream of shock wave 
rocket nozzle throat 


2. Standard Results of Heat and Mass Transfer Theory 
2.1. THE STANDARD PROBLEM 


Figure 1 illustrates the essential problem of mass transfer theory. Material is 
transferred, at a rate mm” (Ib. mass/ft.? h.), across a boundary separating the fluid 
phase under consideration from a neighbouring phase. 


The question which must be solved is: How is m” related to the state of the 
fluid in the main stream (state G), the state of the fluid prevailing adjacent to the 
boundary (state S), and the state of the transferred substance (state T)? 
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FLUID STREAM 


FLUID PHASE 


state T~ 
TRANSFERRED 
SUBSTANCE 


FiGurE |. Notation used in describing mass transfer across a phase boundary. 


State T is defined as having the net composition (irrespective of the state of 
chemical aggregation) of the material crossing the control surface L and an enthalpy 
such that q”, (B.t.u./ft.? h.) is the only heat flux crossing the fictitious control 
volume formed by the L and T surfaces. 


2.2. THE “OHM’S LAW” OF MASS TRANSFER 


It may be shown that, provided that the flow is turbulent or that a sufficient 
number of diffusivities (material, and perhaps thermal) are equal at every point in 
the flow, the answer to the question just formulated may be put in the form 


where B is the dimensionless driving force for mass transfer, and g is the surface 
conductance for mass transfer. 


2.3. THE DrivING FoRcE 


The driving force B is dependent only on the fluid states just referred to; it may 
be written as 


(2) 


where P is any conserved property, and suffixes G, S, and T denote the fluid state 
for which P is to be evaluated. 


Two conserved properties which will be particularly important in the present 
paper are: the specific enthalpy, A (B.t.u./Ib. mass) of the fluid mixture; and the 
mass fraction, f, in the mixture of material from one of the two streams contributing 
to the mixing field. 
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(4) 


2.4. THE SURFACE CONDUCTANCE 

The conductance, g, is mainly determined by fluid-mechanical quantities. For 
the purposes of the present paper, we may state that ¢ is related to the Stanton 
number of heat transfer theory, (St), by 


where G is the mass velocity (velocity x density) of the main stream (lb. mass/ft.? h.). 
The Stanton number is dimensionless. For forced convection it is a function of 
the Reynolds and Prandtl numbers, the latter being the ratio of the kinematic 


viscosity to the thermal diffusivity of the fluid. It is also, as will be shown, somewhat 
dependent on the value of the driving force, B. Thus 


wherein the function ¢(. .. .) depends upon the geometry of the boundary of 
the stream. 


2.5. HEAT TRANSFER WITHOUT MASS TRANSFER 


If the boundary between the phases is impervious to mass, the heat transfer 
rate across the L control surface reduces to 


where g* is the value of g valid for B=0. 


Since, in addition, f has the same value in the G- and S-states (same composi- 
tion, though possibly differing chemical aggregation), the values of h, and hs are 
unique functions of the temperatures and pressures if thermodynamic equilibrium 
can be presumed to prevail. 


2.6. THEORETICAL PREDICTION OF THE SURFACE CONDUCTANCE 


When the S- and 7-states are uniform over the whole boundary, g can be 
obtained from the solution of the partial differential equation of conservation : — 


G.(Vb)- V{y(Vb)}=0 . ‘ ‘ 
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_ fo-fs 
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The boundary conditions are as follows : — 


At the surface, b=b;=0 (9) 
0b 

and {3} .. | 

In the mainstream, b=b,=B. abn 


Clearly, from equations (1) and (10), g is obtained from the gradient of the potential 
b at the boundary, through 


In equation (8): — 


G is the local vectorial mass velocity (Ib. mass/ft.? h.); it must be 
obtained by solution of the equations of motion and continuity; 


and y (Ib. mass/ft. h.) is equal to diffusion coefficient multiplied by density, 
or to thermal conductivity divided by specific heat at constant 
pressure, whichever is appropriate. 


Solutions of the system of equations and boundary conditions are introduced 
below (Figs. 14, 15 and equations (21), (35)). 


2.7. REPRESENTATION OF FLUID STATES ON ENTHALPY-COMPOSITION DIAGRAMS 

With the restrictions which have been imposed (fully turbulent flow or equal 
diffusivities; uniformity of S- and T-states), a simple representation of the fluid 
states, and of the driving force, becomes permissible, employing a diagram 
commonly used by chemical engineers, the enthalpy-composition chart. 


Figure 2. Enthalpy-composition (h-f) diagram: representation of G-, S- and T-states. 
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The co-ordinates of such a diagram are h and f (Fig. 2). Every equilibrium 
fluid state is represented by a point on the diagram. In particular the important 
mainstream (G), surface (S), and transferred-substance (7) states can be so 
represented. It is a consequence of the theory of mass transfer that the three 
corresponding points lie on a straight line. 


Further, as implied by equations (3) and (4), the driving force is represented by 
a length ratio, namely, 


(13) 


S normally lies on the line representing equilibrium between the fluid and the 
substance on the other side of the phase boundary. The other lines drawn on Fig. 2 
represent isotherms (lines of constant temperature). 


Another important implication of the theory is that all the states within the 
boundary layer are to be found along the line GS. If the fluid at any point is in 
thermodynamic equilibrium, the corresponding temperature can be read from the 
appropriate isotherm. 


2.8. JUSTIFICATION OF THE STATEMENTS MADE 


The theory justifying the foregoing statements has been developed by several 
authors. No complete exposition of the theory in the present terms is currently 
available, but partial accounts have been published“:” and a more complete 
theoretical derivation is in course of publication®; Ref. 3 contains a literature 
survey. 


SURFACE OF 
MISSILE NOSE 


SHOCK WAVE 


Ficure 3. Sketch of shock wave and flow directions adjacent to nose of a hypersonic missile. 
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3. The Laminar Axi-Symmetrical Stagnation Point 
3.1. DESCRIPTION OF THE PROBLEM 


The first of the two geometrical systems to be discussed in detail is the forward 
stagnation point of an axi-symmetrical body. In particular we shall have in mind 
the cooling problem which arises when the flight speed of the body is hypersonic. 


Figure 3 illustrates the flow pattern in the neighbourhood of the nose. A shock 
wave is present a short distance ahead of the solid surface; upstream of this wave the 
relative gas velocity is axial, with value u_.. downstream of the wave, the velocity is 
directed radially outward from the axis of symmetry and has the value ug. The 
quantity u, increases linearly with distance x from the stagnation point, in 
accordance with the equation (according to Li and Geiger“) 


where R,=radius of nose (ft.) 
P-« = density of gas upstream of shock (Ib. mass/ft.*) 
Po= density of gas downstream of shock (Ib. mass/ft.*) 
Pol —> (Cp + — Cv) for very strong shock waves 


Cy, Cy= Mean specific heats at constant pressure and volume respectively 
(B.t.u./Ib. mass °F). 


The right-hand side of equation (14) is ———- equal to 0:55u_,./R. in 
practical circumstances. 


Friction at the solid surface causes the establishment of a thin boundary layer 
of retarded gas adjacent to the surface. This boundary layer can usually be taken 
as laminar, at least in the immediate vicinity of the stagnation point. 


The problem of how to calculate the temperature of the surface of the missile 
nose is now considered. First, it will be supposed that the surface is cooled simply 
by radiation, and by conduction to the interior of the missile; then the use of a 
transpiration-coolant will be considered; finally we consider the process known as 
“ablation.” The second and third cases involve mass transfer, but the first does not. 


3.2. HEAT TRANSFER IN THE ABSENCE OF MASS TRANSFER 
3.2.1. The Surface Temperature 


The surface temperature may be determined by solution of equation (7), the 
left-hand side being given by 


where «=surface emissivity 
o = Stefan’s constant (=0-1713 x 10-* B.t.u./[ft.? h. (°F abs.)*)) 
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T;= absolute temperature of surface (°F abs.) 
q”.=heat flux by conduction to interior of missiie (B.t.u. /ft.? h.). 


3.2.2. The Surface Conductance g* 

To evaluate the right-hand side of equation (7), both g* and (hg—hs) are 
needed. The first of these may be obtained from the solutions of equation (8) for 
isothermal wedges in laminar flow, first obtained by Eckert for a Prandtl number 
of 0:7. The result is, by use of the Mangler transformation which relates two- 
dimensional to axi-symmetrical flows, 


where » and p are respectively “the” viscosity and density of the fluid in the 
boundary layer. 


A more recent calculation by Fay and Riddell” has taken account of the 
variation of viscosity and density with temperature. Using values of the properties 
valid for air in a range of practical interest, they derived a modified form of 
equation (16), namely, 


where suffixes § and G denote respectively the gas in contact with the surface and 
the gas between the boundary layer and the shock wave. 


Clearly, combination of equations (14) and (17), together with insertion of the 
appropriate fluid properties, permits g* to be evaluated in any particular case. 
Incidentally, it is evidently possible to re-write equation (17) as 


which brings out the relatively weak influence of the property variations if 
G-properties are used on the left-hand side; for, since » 0c T°’ and pocT~ for 
common gases, the bracket on the right-hand side of (18) is equal to (T,/Ts)"™. 


3.2.3. The Enthalpy Difference (hg —hs) 


The “air” between the shock wave and the boundary layer is usually highly 
dissociated. This, however, can be entirely ignored in calculating he; for an 


114 The Aeronautical Quarterly 


(17) 


HEAT AND MASS TRANSFER 


Vi Ficure 4. Enthalpy- 
io 


temperature diagram 
for air (Ref. 8). 


100 


h=0-24T 


6,000 72000 €,000 9000 10,000 


application of the Steady-Flow Energy equation to a control volume enclosing a 
section of the shock wave near the axis leads to 


- . @& 


where A_,,=enthalpy of the air far upstream of the shock wave (B.t.u./Ib. mass) 


g,=the constant in Newton’s Second Law of Motion, =4:1 x 10° 
(ft. lb. mass/Ib. force h.’), 


and J= mechanical equivalent of heat, = 778 (ft. lb. force/B.t.u.). 
The enthalpy difference required in equation (7) thus becomes 


. . . 


Here the difference h_,,— hs can be determined by reference to tables or graphs of 
the enthalpy of air versus temperature and pressure. Fig. 4 indicates the nature of 
the variation; the strong effect of pressure indicates the presence of dissociation. 
Usually, however, both h_. and hs correspond to temperatures which are sufficiently 
low for dissociation to be negligible; then the Gas Tables of Keenan and Kaye"? 
can be used. 
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3.2.4. The Heat Flux 


The final expression for the heat flux is now obtained by combining equations 
(14), (18) and (20). This can be solved simultaneously with equation (15) if 
information is available concerning the conductive heat flux. The matter is not 
pursued further here; however, it is worth noting that, when the kinetic energy term 
outweighs the other terms appearing in equation (20), the heat flux to the surface is 
proportional to u_,.°/*, neglecting secondary effects from variations of uc, pc, and 
so on. 


It should be mentioned that Fay and Riddell also considered the case in 
which the diffusivities of heat and matter were not equal (Lewis number not equal 
to unity). This case is not discussed here, since it leads to reaction-kinetic con- 
siderations; it is intended to restrict this paper entirely to processes which can be 
computed from physical considerations alone. 


3.3. TRANSPIRATION COOLING 
3.3.1. General Considerations 


It is often impractical to maintain the surface temperature of a missile nose at 
tolerable temperatures by radiation and internal-conduction cooling alone. An 
alternative possibility is to use transpiration cooling, e.g. to force hydrogen through 
the nose surface, which is provided with small pores. Fig. 5 illustrates this. 


As will be seen, it is desirable that the coolant should have a high specific heat; 
hydrogen gas is therefore suitable. On the other hand, hydrogen is liable to burn 
at the high temperatures which are in question; a flame surface is therefore drawn 
in Fig. 5. 


G-STATE 
_COOLANT STREAM 
S-STATE 


T-STATE 


Ficure 5. Transpiration cooling of a missile nose. 
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Ficure 6. Enthalpy-composition 
- diagram for H,-air at low pres- 
~N é sures, used for determination of 


« B (Ref. 10, corrected). 


1000 °F 


The theoretical problem which must be solved by the designer is: How much 
coolant must be forced through unit area of the surface in unit time in order to 
maintain a prescribed surface temperature? This problem can be solved by the use 
of standard techniques of mass transfer theory, as will now be demonstrated. 


3.3.2. Determination of the Driving Force 

The required coolant flow rate can be regarded as a mass transfer rate, m”, to 
be evaluated from equation (1). For this a conductance g and a driving force B 
are needed. 


The driving force B will be determined first, using the method of Ref. 10. We 
need an enthalpy-composition (h-f) diagram for hydrogen-air mixtures. This is 
shown in Fig. 6; the lines on it are isotherms. 


Making use of the results of Section 2.7, we note that the point G (representing 
the G-state) must lie on the vertical line f=0 (pure air), at a height given by hz 
calculated from equation (19). The point T (representing the 7-state, i.e. pure 
hydrogen at the temperature of the coolant reservoir), lies on the vertical line, f=1, 
and the appropriate isotherm. The point S is therefore located as the intersection 
of the line G7 and the prescribed ft, isotherm.* 


° that thermodynamic at the It need not prevail anywhere 
else, however. 
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Ficure 7. The func- c 4 
tions E(B) and F(B) 
for use with equation o3 7 40 


(21) for laminar 
boundary layers with 
(Pr)=0°7; also, 
+2) for use with of /F+2)- 
equation (22) for the 
axi-symmetrical stag- 


Having located G, T, and S in this elementary fashion, equation (13) is invoked 
for the evaluation of B. It is equal to 1-725 in the example shown, in which 
ts=2,000°F abs., h, = 8,000 B.t.u./lb. mass (corresponding to u_,, ~ 20,000 ft. /sec.), 
and the temperature of the coolant reservoir is 500°F abs. If surface radiation is to 
be accounted for, hy is reduced by the magnitude eoT,*/m”. 


3.3.3. Determination of g 


The next step is the evaluation of the surface conductance g. Relatively few 
exact solutions of the system of equations (8), (9), (10) and (11) are available. For 
laminar boundary layers, the most compact approximate calculation formula is that 
of Ref. 11, which may be written in the form valid for a Prandtl number of 0-7 as 


1/2 


x 
| upuc! uc?-*R*dx | (21) 
where E and F are the functions of the driving force B illustrated in Fig. 7. 
Since, for the axi-symmetrical stagnation point, R, the distance from the axis, 
is equal to x, while u, is proportional to x, equation (21) reduces to 


which differs from equation (16) only in the substitution of the quantity E /(F + 2) 
for the constant 0-957. Fig. 7 also contains E/ (F +2) as a function of B.* 


*Note: At B=0, EV (F +2) reduces to 0-923 and not to 0-957, the exact value; this difference is 
a consequence of the approximate nature of equation (21). For the same reason, the curve 
of Fig. 8 is slightly below the more exact one to be found in Ref. 10. 
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3.3.4. Determination of the Coolant Rate 


With g and B calculated, the desired coolant rate m” is obtained from 
equation (1). To make clear the effect of the driving force on sn”, equations (1) and 
(22) have been combined, yielding the curve shown in Fig. 8. It is evident that the 
coolant rate does not increase linearly with the driving force. 


The curve of Fig. 8, together with an enthalpy-composition diagram, such as 
Fig. 6, and evaluation of equation (14), permits the solution of the problem. Only 
standard techniques have been used: none have had to be invented specially for 
this problem. 


The curve of Fig. 8 can also be used, with little error, when the fluid properties 
» and p are inserted at their values just outside the boundary layer, as indicated by 
the suffix G. This is confirmed by the data of Fay and Riddell” and by the more 
recent calculations of Howe and Mersman®*’. Points derived from Ref. 23 are 
plotted in Fig. 8; it can be seen that they lie somewhat above the curve based on 
equation (22), which is, of course, approximate. 


3.4. THe USE oF A GRAPHITE “HEAT-SHIELD” 


3.4.1. General Considerations 


The provision of a porous nose and a correctly metered supply of transpiration 
coolant obviously involves the solution of difficult engineering problems. There is 


therefore much to be said for using a solid nose material which produces similar 
effects by surface vaporisation or melting. This process has come to be called 
“ablation”; it is, of course, a particular sort of mass transfer. 


A possible material is graphite. When sufficiently hot, this material burns 
away when exposed to a stream of air. We will examine how the rate of burning 
can be calculated. As before, two quantities have to be evaluated, the surface 
conductance g and the driving force B. There is nothing new to be said about the 
first of these steps: g can be determined from equation (22), as before. We will 
therefore concentrate on the calculation of B. 


3.4.2. Enthalpy-Composition Diagram 


h-f diagrams for carbon-air mixtures have been published in Ref. 12. One 
valid for a pressure of one atmosphere is shown in Fig. 9. The diagram shows 
isotherms as before; but there is also a new feature: the phase boundary or 
“saturation line”, along which are to be found gas states in thermodynamic equi- 
librium with solid carbon. 


The state-point of the main stream, G, is found as before; it lies on the pure-air 
(f =0) line and has an enthalpy determined from equation (19). A, is 8,000 B.t.u./ 
Ib. mass in Fig. 9, the same value as was used in the previous example. 


The state-point of the transferred substance, T, lies on the pure-carbon (f= 1) 
line; its vertical position is determined by the initial temperature of the graphite 
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Fioure 8. Curve showing mass transfer rate at axi-symmetrical stagnation point. (Pr)=0-7. 
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Ficure 9. Enthalpy-composition diagram for carbon-air mixture at 1 atm. pressure (Ref. 12). 
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mass, 500°F abs. in the example illustrated. If surface radiation were to be con- 
sidered, however, T would be lower by an amount corresponding to the enthalpy 
difference :oT;*/m’”. 


__ The state-point of the gas in contact with the surface, S, must lie on the line 
GT; but it must also lie on the “saturation line” if thermodynamic equilibrium 
prevails at the surface. The location of S is therefore easily found. 


Measurement of the ratio GS/ST now yields the value of the driving force B 
by way of equation (13). B is equal to 0:39 in the example shown. It should be 
remarked that S lies on the sloping part of the saturation line: approximately half 
the carbon is transferred in the form of vapour in the example calculated. The 
surface temperature is about 7,800°F abs. 


3.4.3. Determination of the “Ablation” Rate m” 


Knowledge of the driving force, together with the curve of Fig. 8 and an 
evaluation of equation (14), enables the mass transfer rate m” to be determined. 
From this it is possible to calculate how thick the “heat-shield” must be. 
Evidently, since B for graphite is less than a third of that just calculated for hydro- 
gen, the use of a graphite “heat-shield” should lead to a reduction in missile weight. 


3.4.4. “Heat Stored” in the Shield 


The mass transfer rate just calculated is that valid for the steady state. At 
first the rate will be lower, since the graphite is cold. The transient heating problem 
can be calculated by conventional techniques. 


It is interesting to calculate the “heat stored” within the shield in the steady 
state, because, after the end of the high-velocity flight, this may “soak” into the 
body of the missile. It may be easily shown that this quantity is kp (ts—fo)/m”, 
where k is the thermal conductivity of the graphite, p its density, and f: and f» are 
the surface and initial graphite temperatures respectively. Clearly it is desirable 
to use a graphite of low conductivity and density. Materials which form a porous 
carbon-containing structure on heating are advantageous in this regard. 


3.5. “MELTING ABLATION” 


Consideration has been given by various authors (e.g. Ref. 13) to the use of 
“heat-shield” materials which melt on heating. Although such materials are of 
minor practical importance, it may be worth examining how the use of an enthalpy- 
composition diagram can illuminate the utility of a particular material. 


3.5.1. Enthalpy-Composition Diagram 


Figure 10 represents an hA-f diagram which might be appropriate to such a 
material mixed with air at a fixed pressure. The phases present where the state- 
point of the mixture lies in a given area are marked on the diagram; a few 
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representative isotherms are also shown. The straight lines separating the (gas + 
liquid + solid) region from its neighbours both represent the isotherm corresponding 
to the melting temperature. 


The point G representing the air stream state appears on the f=0 line as before. 
The point O, representing the initial state of the bulk of the shield material, lies on 
the f=1 line at the end of the appropriate isotherm. 


3.5.2. Determining the Mass-Transfer Driving Force 


We shall consider the case in which the film of liquid separating the gas stream 
from the solid surface is so thin as to be at nearly uniform temperature. This 
occurs when the liquid viscosity is very small or its thermal conductivity very high. 
An immediate consequence of the assumption is that the state-point § must lie in 
the position shown on the diagram; for the gas in contact with the surface is in 
equilibrium both with liquid and solid shield material. 


With S located, the state of the transferred substance is also found; for T lies 
on GS produced, and also on the line f=1 (assuming that the melted material is 
unoxidised). The mass transfer driving force, B, is equal to the ratio GS/ST, as 
before. Knowledge of B then leads to that of m”, the rate of mass transfer into the 
gas phase. 


40) 
GAS FSQLID 
( = 
100% PURE AIR f ngs | 
100 ABLATING 
SUBSTANCE JT 


Ficure 10. Enthalpy-composition diagram for a hypothetical substance exhibiting 
“melting ablation”. 
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Figure 11. Control volume for analysing mass transfer and melting at a solid surface. 


3.5.3. Calculation of Total Ablation Rate 


Of major interest to the designer is the total rate at which shield material 
disappears, m”’,.:; this comprises both mm” and the rate at which material melts away. 


The state-point of the melted material is L on Fig. 10. An application of the 
Steady-Flow Energy Equation to the adiabatic* control volume of Fig. 11 therefore 
evidently results in 


(23) 


where the lengths appearing on the right-hand side obviously refer to enthalpy 
differences. 


A short and obvious train of argument thereupon shows that if a line is drawn 
through G, parallel to SL and to cut OL at U, then 


(24) 


3.5.4. Practical Consequences 


Two conclusions can be drawn from equation (24). The first is that the 
melting rate increases with LU, and so with the magnitude of the gas stream 
enthalpy, and decreases with OL, i.e. with the heat required to melt unit mass of 
the material. The second conclusion is that, if a new material is to be assessed for 
its suitability to be used in a heat-shield, it is very important to have the necessary 
thermodynamic data. It should be clear from Fig. 10 that these include knowledge 


*An obvious modification needs to be made if surface radiation is important. 
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of the gas-liquid equilibrium when mixed with air, as well as knowledge of the 
solid specific heat and of the latent heat of fusion. It is hard, in the author’s 
opinion, to find a better way of displaying such data than representation on an 


enthalpy-composition diagram. 


3.5.5. Related Processes 


We shall not consider the somewhat more complex case in which a finite 
temperature difference exists across the liquid film, beyond remarking that this 
phenomenon leads to a reduction in the melting rate, and that similar problems 
arise in, for example, steam condensers; calculation techniques have been developed. 
The melting ablation process has features in common with the cutting (by burning 
and melting) of steel plates by oxygen jets, and this process can be treated from 
the point of view of mass transfer theory; agreement between theory and experiment 
is 


Sometimes the mass transfer into the gas phase proceeds by way of a tem- 
perature-dependent pyrolysis reaction; then chemical-kinetic data are needed. 
When these data are available, it is easy to devise appropriate modifications to the 
calculation procedure for the ablation rate. Similar problems arise in the theory 
of solid-propellant combustion. 


G, 
Lok =( = D 


KEROSINE 


Fioure 12. Sketch of rocket motor burning liquid oxygen and kerosine. 


4. Heat and Mass Transfer at the Throat of a Rocket Nozzle 
4.1. DESCRIPTION OF THE PROBLEM 

Attention is now turned to the other end of the missile, the rocket motor. 
For concreteness, we shall consider a regeneratively cooled LOX*-kerosine motor 
of the type sketched in Fig. 12. A major problem for the motor designer is the 
cooling of the wall of the nozzle, particularly in the neighbourhood of the throat. 
We shall consider how the theory of mass transfer can assist in the prediction of 
heat transfer rates. 


*LOX=liquid oxygen. 
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Four cases are considered. Firstly we shall inquire how to calculate the heat 
transfer rate, assuming that the bare metal wall is exposed to the hot gases. We 
then consider film-cooling, i.e. that method of cooling in which liquid fuel is sprayed 
on to the wall at the injector end and forms a continuous film at least over the 
walls of the combustion chamber and of the upstream portion of the nozzle. 


In practice, a film of another kind often forms on the nozzle wall, a film of 
carbon. This has the valuable effect of acting as a heat insulator. We shall 
attempt to understand the laws governing its thickness. Finally, the rate of burn- 
out of an uncooled graphite nozzle throat will be considered. 


In all these cases the problem presents the same two aspects which have been 
encountered before: the determination of a driving force and the determination of a 
conductance. For the former, enthalpy-composition diagrams will be used; for the 
latter, we shall make use of boundary layer theory. This time however the 
boundary layer is turbulent. 


4.2. Heat TRANSFER TO A METAL NOZZLE THROAT 
4.2.1. Approximate Theory 


An upper limit to the heat flux can be calculated extremely simply. Equation 
(7) can be written as 


where (St)* is the Stanton number for zero mass transfer, which usually has a value 
around 0-002, G is the mass velocity in the nozzle, and h,,°, the stagnation enthalpy 
of the gas stream, has been substituted for h,, because high Mach numbers are in 
question. 


Clearly G has its maximum value where the nozzle has its minimum area, i.e. at 
the throat. We will concentrate on this region, and call the mass velocity there G,. 
Now G, is related to the characteristic velocity c* of the propellant mixture, 
measured in ft./sec., by 


G.=2-45 x 10° Ib, mass/ft.? h., 


where p, is the combustion chamber pressure in atmospheres and c* is a tabulated 
property of the propellant mixture; for LOX-kerosine c* has the value®” 
5,540 ft./sec. 


Now h,° —hs cannot exceed Ah,,;:, the enthalpy difference between the gas in 
the combustion chamber and that when the gases have expanded through an infinite 
pressure ratio. But Ahn: is related to c* by 


& 9:25 10-° B.t.u./Ib. mass. 
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In equation (27) c* is measured in the same units as before, and the maximum value 
of the thrust coefficient has been taken as equal to 2°15. 


Inserting (St)* =0-002, we obtain a simple formula for the heat flux at the 
throat, namely, 


mx ™~45-3p,c* Btu/ft?h . (28) 


A LOX-kerosine motor operating at 20 atmospheres combustion-chamber pressure 
therefore has an upper limit to the throat heat flux given by 


max = 45°73 x 20 x 5,540 B.t.u./ft.2 h.=5-02 x 10° B.t.u./ft.? h. 
B.t.u./in.? sec. 


In practice heat fluxes of about one half this magnitude are encountered, as might 
be expected. 


In addition to giving an upper limit to the heat flux at the throat, equation (28) 
also indicates the magnitude of the effect of a change in combustion-chamber 
pressure or a change of propellant: the heat flux is proportional to the pressure and 
to c* which measures the “energy” of the propellant. It is worth noting that the 
whole range of c* for common propellants is not large: of 23 bi-propellant 
combinations tabulated by Sutton”, the highest c* value is 7,930 ft./sec. (for 
LOX-liquid-hydrogen), and the lowest is 4,940 ft./sec. (HTP-C-Stoff). 


4.2.2. More Precise Determination of Enthalpy Difference . 


Returning now to equation (25), we consider how the term (h,,°—hs) may be 
evaluated. Once again, an enthalpy-composition diagram is used. 


Figure 13 shows an h-f diagram for LOX-kerosine mixtures at 20 atmospheres 
pressure. It will be seen that there are now more than one saturation lines, since 
both fuel and oxidant exhibit liquid phases. The state-points F and O represent 
respectively the kerosine and liquid oxygen at entry to the chamber. 


If the flow and chemical-reaction processes in the combustion chamber are 
regarded as comprising a steady-flow adiabatic mixing process, the laws of thermo- 
dynamics imply that the stagnation condition of the combustion gases is represented 
on the diagram by a point G° lying on the line FO; moreover its location is such that 
the length ratio FG’/G°O is equal to the oxidant-fuel ratio at which the motor 
operates. In Fig. 13, G° has been drawn to lie on the line f=0-3, which 
corresponds roughly to practical conditions. 


The point S, representing the gas in contact with the metal surface, must have 
the same f value as G’; further, it must lie on the isotherm corresponding to the 
metal temperature ts, which will be supposed known. These considerations locate 
S, and so enable h,°—h, to be evaluated as the vertical distance SG. In the 
example shown in Fig. 13, ts=1,000°R and h,° — hs =2,625 B.t.u./Ib. mass.* 


*Note: Strictly, we should use the ¢, isotherm on a diagram corresponding to the pressure at 
the section of the nozzle in question. Fortunately the influence of pressure is very small at the 
wall temperature so that this refinement is unnecessary. 
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Figure 13. Enthalpy-composition diagram for LOX-kerosine mixtures at 20 atm. pressure. 


-5000 


It should be noted that no inquiry has had to be made into the degree of 
dissociation present within the combustion chamber nor into the question of what 
mean specific heats are appropriate. Provided that the basic assumptions hold 
(equal diffusivities), all such difficulties are circumvented. 


4.2.3. More Precise Determination of the Stanton Number 


Methods for the determination of the Stanton number are not very advanced. 
Sometimes the pipe-flow formula is used“, namely, 


(St)* =0-023 : (29) 
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In equation (29), 
D=nozzle diameter at section in question (ft.) 
=gas viscosity (Ib. mass/ft. h.) 
(Pr)=gas Prandtl number. 


However this practice is undoubtedly too crude; for the accelerating flow results in 
quite a thin boundary layer. 


A more sophisticated practice is to use the theory of the velocity boundary 
layer to calculate the distribution of drag coefficient c; by the way of the von 
K4rm4n momentum equation. Thereafter, (St)* is found by the application of 
some form of the Reynolds Analogy 


However, this is also inaccurate and is known to overestimate (St) in accelerating 
flows. The continued use of this practice is surprising, since it has long been 
recognised in connection with laminar heat transfer that, if only one boundary layer 
is to be considered, it is better to deal directly with the thermal boundary layer. 


The latter technique has been proposed for turbulent boundary layers by 
Ambrok“®. Adapting his procedure to the requirements of compressible flow in an 
axi-symmetrical duct of diameter D, we obtain, for a fluid Prandtl number of 0-7, 


where x is the distance along the nozzle wall from the injector and «, the gas 
viscosity, is to be evaluated at some value intermediate y and ps. 


Equation (31) is easily evaluated, since D and G are known as functions of 
distance. For rockets, for which 
(32) 


where A is the cross-sectional area of the nozzle section, and suffix “‘t’’ refers to the 
throat, equation (31) simplifies to 


(D,/D)® (33) 


2/D, 
0-75 py 


This expression is almost as easy to evaluate as the pipe-flow formula (29) 
and rests on a much more secure foundation. The quadrature is of course entirely 
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concerned with the motor geometry. Since » is usually around 0-1225 Ib. mass/ft. h. 
under rocket motor conditions, equation (26) permits the Reynolds number term 
to be put in the simplified form 


where p, is in atmospheres, D, is in ft., and c* is in ft./sec. 


4.2.4. Determination of the Heat Flux 


The convective heat flux can now be calculated by combination of equations 
(7), (26), (33) and (34). It is easy to see that the maximum heat flux occurs very 
near the throat, because of the D,/D)''* term in the numerator of equation (33). 
As an example, we calculate the throat heat flux of a LOX-kerosine motor with a 
1 ft. throat diameter operating at 20 atmospheres pressure with the G° and S§ points 
of Fig. 13; we assume that the quantity 


z/D, 
D. 0-75 dx 


in equation (33) is unity, as it often is. There results 
= 6°95 B.t.u./in.? sec. 


This heat flux is lower than the maximum value calculated earlier, but it is 
higher than values commonly encountered in practice. This is largely because the 
gases near the wall are more fuel-rich than the average; but another possible reason 
for the discrepancy, namely the carbon film, is considered later. 


4.3. WAaLL-CooLING BY MEANS OF A FUEL FILM 
4.3.1. Statement of the Problem 


It is possible, in principle, to dispense with regenerative cooling if a continuous 
film of liquid fuel can be maintained over the combustion chamber and nozzle 
walls. The liquid might be sprayed on at the injector end, or it might be caused 
to seep through pores in the walls. The question of how much fuel vaporises per 
unit time per unit area of film is now examined, for this question must be answered 
when the feasibility of the method is being considered. 


This is a straightforward mass-transfer problem. To solve it, it is necessary 
to determine a driving force and a surface conductance. 


4.3.2. The Driving Force 


It is easily shown''”’ that the driving force is determined by locating a state- 
point S; on the fuel-saturation line of the enthalpy-composition diagram such that 
the line G°S; coincides with a mixed-phase isotherm. Fig. 13 illustrates this. With 
the transferred-substance state-point T lying on the line f=1, the driving force B 
is found from equation (13); in the example shown, B = 10. 
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4.3.3. The Surface Conductance 


To calculate g, we use a modification of equation (31), namely, 


x 


The function of B appearing on the right-hand side is only approximate, but it fits 
the more complex function calculated by Rubesin''” fairly well, and the experi- 
mental data of Mickley and Davis''*’ equally well. 


4.3.4. Determination of the Vaporisation Rate 


With B and g calculated, it merely remains to substitute in the “Ohm’s Law” 
equation (1), and so to deduce the distribution of local vaporisation rate, m”’. A 
quadrature over the whole surface then leads to the total mass flow rate of fuel 
which must be supplied to the wall. 


Many other questions relating to film-cooling need to be answered by the 
designer, for example: —‘Will the liquid film be ruptured by the aerodynamic 
forces acting on it?” “Does the unevenness of fuel/oxidant ratio across the nozzle 
lead to an appreciable performance reduction of the whole motor?” “What happens 
when the combustion-chamber pressure exceeds that of the critical point of the 
fuel?” These questions are not discussed here. 


4.4. THe INFLUENCE OF CARBON DEPOSITION ON THE WALL 


As mentioned already, the actual rates of heat transfer encountered in LOX- 
kerosine motors are appreciably lower than those calculated by methods such as 
that of Section 4.2. The reason is believed to be that a film of carbon builds up on 
the wall and acts as a thermal insulator. The conditions under which this is likely 
to occur, and the magnitude of its effect are now discussed. 


4.4.1. The Carbon-Saturation Line 


The enthalpy-composition diagram for LOX-kerosine (Fig. 13) shows a broken 
line marked “carbon saturation”. This connects the states at which solid carbon 
is just about to be precipitated from the mixture: to the left of the line, there is no 
free carbon; to the right of the line, soot particles will form if equilibrium is 
attained. 


The state-point S, representing the gas in contact with the metal in our example, 
lies below the carbon-saturation line. This means that, if metal at this temperature 
were exposed to the gas, carbon would condense on the metal; the process is akin 
to frost formation on refrigerant coils. 
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Without analysing the transient process of carbon build-up, we can see that a 
steady state can be reached when the gas in contact with the carbon film has the 
state-point shown as S. in Fig. 12. The corresponding temperature is higher than 
that of S; the thickness of the carbon film must be such as to produce the 
corresponding temperature difference by way of its “resistance” to heat flow: this 
temperature difference is nearly 1,000°F in the example shown. 


4.4.2. The Heat Flux 


___ It is clear that the heat flux to the wall has been decreased in the ratio 
G°S./G°S, about 0-8 in the case shown. Now it is common practice to design the 
injector so that the gas in the immediate neighbourhood of the motor wall has a 
higher value of f than the average; thus the state point G’ which is relevant to heat 
transfer may be shifted to the right. Because of the sieepness of the carbon- 
saturation line, it is evident that quite a small shift can produce a large decrease in 
the enthalpy difference causing heat transfer, represented by G’S.; the carbon film 
increases in importance the farther G moves to the right. The thickness of the 
carbon film is easily calculated from the heat flux if the thermal conductivity of 
carbon is known. 


In the absence of quantitative experimental confirmation, including measure- 
ments of carbon-surface temperatures, film thickness, and so on, the present 
explanation must be regarded as speculative. However, it is entirely in accordance 


with the general theory of mass transfer and with experience of other types of 
equipment, for example condensers. 


The carbon film will be subjected to rather great shear stresses as a result of the 
high gas speed; it is therefore likely to be torn off periodically. This should be 
remembered before the existence of the film is relied on too strongly in design. 


4.5. THE BURN-OUT OF A CARBON THROAT 


Solid-propellant rockets cannot be regeneratively cooled; the nozzle material 
must therefore be capable of withstanding high temperatures. Graphite may be 
considered for this purpose. As our final problem, we shall discuss how to 
calculate the rate at which the carbon is likely to burn away during the course of a 
firing. Attention will be concentrated on the throat, since it is there that diameter 
changes can have their greatest effect on the performance of the motor as a whole. 


4.5.1. Procedure for Calculating the Burn-Out Rate 


The steps in the calculation are precisely the same as before. Firstly, the 
problem is recognised as one of mass transfer. Secondly, an enthalpy-composition 
chart is constructed for mixtures of carbon with the combustion products of the 
propellant in question; strictly speaking, several such charts are needed, valid for a 
range of pressures. A graphical construction on the A-f diagram then leads to an 
evaluation of the appropriate driving force for mass transfer. 
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The next step is to determine the mass transfer conductance g by evaluating 
equation (35) for the configuration, pressure and propellant properties in question. 
Finally, the burn-out rate rm” follows by evaluation of equation (1). 


4.5.2. Analytical Determination of the Driving Force 


Complete enthalpy-composition charts of the type required are seldom 
available. In their absence, it is usually sufficiently accurate to calculate B from 
the composition of the combustion products by way of a formula derived in Ref. 2, 
namely, 


where mo,,¢, Mco_c, and so on, are the mass fractions, in the main gas stream, of the 
gases referred to. 


We shall consider, as an example, the double-base propellant JPN Ballistite, of 
which the combustion products at 5,500°F and 1,500 Ib./in.? are described, 
according to Ref. 19, by 


mo,,¢=0-006; moo,,¢=0°352; mu,o,e=0°175. 
Insertion of these data in equation (36) yields B=0-217. 


4.5.3. Rough Determination of the Burn-Out Rate 
A similar theory to that underlying equation (28) leads to 


ri”, =4-9 x 10° p,B/ck Ib. mass/ft?h, . .  . (37) 


where p, is the combustion-chamber pressure in atmospheres and c* is the 
characteristic velocity of the propellant in ft./sec. 


Inserting the values relevant to JPN Ballistite at p,=100 atm. (c* ~ 4,000 
ft./sec.), we obtain 


m”.=2°66 x 10° Ib. mass/ft. h. 


Since the density of graphite is about 142 Ib. mass/ft.*, it is seen that the throat 
radius will increase at the rate of 0:0625 in./sec. Whether this “erosion” rate is 
tolerable depends upon the initial throat diameter and the total time of firing. 


5. Outstanding Tasks in Mass Transfer Theory 


It should by now be apparent that standard techniques of mass transfer theory 
can be applied to a variety of aeronautical problems, and that the prediction of the 
corresponding heat and mass transfer rates can be made with remarkably little 
calculation. However, lest it should be concluded that the theory therefore requires 
no further development, some of its weaker aspects will now be indicated. 
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Ficure 14. Surface conductance for a flat plate in laminar flow; exact solutions of the boundary 
layer equations (Ref. 20). «/y= Prandtl or Schmidt number. 


5.1. MASS TRANSFER THROUGH LAMINAR BOUNDARY LAYERS 


Equation (21), with Fig. 7, provides a general method of calculation for mass 
transfer in laminar boundary layers; it represents an interpolation formula based on 
exact solutions of the equations for laminar flow on wedges. However, these 
solutions are extremely scanty; almost all existent solutions can be represented on 
two diagrams (Figs. 14 and 15). The first of these relates to the flat plate; it will be 
seen that values for high and low Prandtl number are lacking. The second is valid 
only for a Prandtl number of 0:7; only positive values of B can be represented and 
a quite small range of Euler numbers. 


It is clear that much work remains to be done before the blank areas on the 
charts have been filled out. This is work which aerodynamicists are well fitted to 
undertake. Perhaps they will be additionally encouraged to do this work by the 
fact that the results are of wide application in mechanical and chemical engineering 
as well as in aeronautics. 


5.2. MASS TRANSFER THROUGH TURBULENT BOUNDARY LAYERS 


Equation (35) represents a general method for calculating mass transfer rates 
through turbulent boundary layers with a Prandtl number of unity. It is, however, 
simpler than the facts warrant: at least the constant on the right-hand side should 
be replaced by a function of Mach number and the ratio between the wail 
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Ficure 15. Surface conductance for similar laminar flows with «/y=0-7; exact solutions of 
boundary layer equations (Refs. 20, 21, 22). 


temperature and the stream temperature; and the relative thicknesses of the velocity 
and compositional boundary layers ought to be taken account of. The B-function 
is also no more than an approximation. 


Although equation (35) could be slightly improved by a full incorporation of 
existing knowledge, we cannot yet quantitatively express the influences just referred 
to. This is a field in which research is required urgently; once again, it is the 
aerodynamicists who have most familiarity with the relevant research techniques, 
both theoretical and experimental. 
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5.3. EFFECTS OF UNEQUAL DIFFUSIVITIES AND OF REACTION KINETICS 


It must be remembered that most of the results discussed hold only if certain 
equalities exist between the transport properties of the fluids; only then can 
reaction-kinetic considerations be neglected. The results are therefore all, to a 
greater or lesser extent, approximate. Clearly the extent of their inaccuracy needs 
to be examined. 


The author does not think this matter so urgent as the research tasks of 
Sections 5.1 and 5.2. We already know roughly what are the effects of the 
influences we neglect, and can devise methods for taking full account of them when 
required; however, at present there does not seem to be any prospect of reducing the 
results of the more exact studies to a form simple enough for the designers to use; 
for, if the equations which must be solved become much more complex than those 
in the present paper, few people will have enough time to use them. 


5.4. FinaL REMARKS 


Probably the greatest stimulus to the development of the general heat- and mass- 
transfer techniques would be their intensive employment in the solution of the 
problems which have been discussed. Thermodynamicists need to make a contri- 
bution here: for enthalpy-composition diagrams have not yet been constructed for 
the range of mixtures and pressures which are encountered in practice. In many 
cases the data are available in the form of tables and scattered reports; but they 
need to be collected and expressed graphically. In some cases the attempt to do 
this will show up gaps in the available data; then experimental research will be 
necessary. 
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The Bending and Twisting of a Flat Elastic 
Plate with a Reinforced Circular Hole 


R. HICKS 
(G.E.C. Simon-Carves Atomic Energy Group) 


SumMARY: This paper considers the problem of a reinforced circular hole 
in a bent or twisted flat elastic plate having small deflections. The reinforce- 
ment is chosen so that the ratio of its radial width to the diameter of the hole 
is small; this implies that there is no radial stress variation across sections of 
the reinforcement. Using the classical theory for bending, expressions are 
obtained for the boundary conditions at the junction of the reinforcement and 
plate. The expressions are general in so far as they are also applicable to 
plates under the action of lateral loading. The following particular problems 
are considered :—(a) uniform bending about one or two axes, (b) uniform 
twisting about two axes. For the particular case of uniform bending about 
One axis, data has been obtained for the stress concentration around holes 
for plates and reinforcements having a reasonably wide range of 
practical dimensions. 


1. Introduction 


It has been shown by Goodier"? that the presence of an unreinforced circular 
hole in a bent or twisted plate results in an increase of the maximum moment. For 
example, the maximum moment in a perforated plate uniformly bent about one 
axis is approximately 70 per cent greater than the maximum moment in a similarly 
loaded plate with no hole. 


Gvodier also considered the effect of the relative size of the hole on the moment 
distribution in the plate. He found that when the ratio of the width of the plate to 
the diameter of the hole is reasonably large (say > 6), the change in the moment 
distribution due to the introduction of the hole is local and results in a stress 
concentration which is independent of the relative size of the hole. 


Goland® investigated the influence of several types of inclusions on the stress 
distributions in elastic plates under transverse flexure. In particular he obtained 
solutions for circular inclusions with various degrees of elasticity and for rigid 
inclusions with differing elliptical forms. 


The problem of an isotropic plate bent or twisted by couples at infinity, and 
containing an unreinforced hole of fairly general shape, has been considered by 


Received July 1959. 
May 1960 


R. HICKS 


Holgate®. He obtained numerical values for the stress couples around the edges 
of square and triangular holes with rounded corners for plates under the action of 
several applied load systems. 


Subsequently, for bent or twisted plates which may be regarded as infinitely 
extended in two directions, solutions for a variety of problems have been found by 
different authors“, In all cases, the solutions are for plates with an unreinforced 
hole or a solid elastic or rigid inclusion. Also a theoretical solution for an unrein- 
forced circular hole in an infinite strip has been obtained by Tamate™. 


Previously, the author“® '” considered the problem of a reinforced hole in an 
end-loaded infinite flat plate. In this case, it was assumed that the reinforcement is 
placed on one side of the plate only, so that the other side remains flush. That is, 
the plate is stressed by the simultaneous action of stresses applied at its middle 
plane and by bending stresses induced locally by the asymmetry of the reinforcement. 
It was found that the induced bending produced stress concentrations which were 
considerably greater than the stress concentrations associated with an identical 
symmetrically placed reinforcement. 


In this discussion a solution is obtained for the stress distribution in a bent or 
twisted plate with a reinforced circular hole. The reinforcement is chosen so that 
the ratio of its radial width to the diameter of the hole is small; this implies that 
there is no stress variation across radial sections of the reinforcement, so that its 
bending and twisting resistance is the same as that of an elastic ring. Furthermore, 
it is assumed that the reinforcement has a symmetrical radial cross section so that 
its neutral axis coincides with the middle plane of the bent plate. 


Using the classical theory for bending, expressions are obtained for the 
boundary conditions at the junction of the reinforcement and plate and for the 
moments and shear forces induced in these elements. The expressions are general 
in so far as they are also applicable to plates under the action of a lateral load 
system, provided that in the region near the hole the deflection surface of a similarly 
loaded unperforated plate can be expressed in polar co-ordinates. 


Solutions for various types of lateral loading are outside the scope of this paper, 
but the following important problems are considered in detail : — 


(a) Uniform bending about one or two axes, 
(b) Uniform twisting about two axes. 


For the particular case of uniform bending about one axis, data has been 
obtained for the stress concentrations around holes for plates and reinforcements 
having a reasonably wide range of practical dimensions. Using the Principle of 
Superposition, this data can also be used to find stress concentrations for other types 
of bending and twisting moment distributions. 


Finally, it is shown that suitable reinforcements can be designed to reduce the 
stress concentration in the plate; however, badly designed reinforcements can result 
in stress concentrations which are considerably higher than those in a similarly 
loaded plate with an unreinforced hole. 
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NOTATION 
r,@ polar co-ordinates 


x,y Cartesian co-ordinates 
w deflection of perforated plate 
w, deflection of unperforated plate 
M.,Me,M. bending and twisting moments in the plate 
M,,M, applied moments bending the plate about the x and y axes 
M, moment applied normally to the edge of a triangular plate 
M,T bending and twisting moments in the reinforcement 
M.,M,y components of M due to change in curvature and radial rotation 
T;,Ts components of T due to twist and tangential slope 
Q,.Qs shear forces in the plate 
V,,Vs shear forces applied to an annular plate 
S shear force in reinforcement 
fz bending stress in reinforcement due to a radial rotation 
a_ radius of hole 
h_ thickness of plate 
t radial width of reinforcement 
d_ depth of reinforcement 
z distance from neutral axis to an arbitrary point in the reinforcement 
Xr twist per unit length of reinforcement 
Xs twist per unit length due to tangential slope of reinforcement 
D flexural rigidity of plate 
1 moment of inertia of reinforcement 
K twisting constant for reinforcement 
y,A4 bending and twisting constants defined in equations (22) 
E Young’s modulus 
Poisson’s ratio 
S length parameter 
C,(r),S.(r) functions of r 


2. Theory 


Figure 1 shows an elemental length of reinforcement subtended by an angle 66 
at the centre of the hole. M and T are the bending and twisting moments acting on 
a radial cross section of the reinforcement and S§ is the accompanying shear force. 
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Ficure 1. Bending moments and shear forces acting on an elemental length of reinforcement. 


This shear force acts in a direction perpendicular to the surface of the bent plate 
and produces the moment S (a — 1/2) 56 which tends to bend the reinforcement about 
a radial line in the plane of the plate. 


(M,), ad0 and (M,), a5@ are the bending and twisting moments transmitted to 
an elemental length of reinforcement by the bent plate, while (Q,),a5@ is the 
corresponding shear force. There are no horizontal shear forces and no forces 
normal to the sides of the element, since strains at the middle plane of the plate and 
reinforcement are negligible when the deflections are small. 


In this analysis the deflection surface of the bent plate is expressed in polar 
co-ordinates having their origin at the centre of the hole. Downward deflections 
are considered positive, so that M, etc. can be written as follows. 


1 dw 1 


2. 
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0 1 dw 1 

D dw 1 
In addition to these relationships, there are the following relationships which give 
the resultant shear on radial and tangential boundaries of the plate : — 


v.=(0-1 5 Me) 


These relationships, which were first obtained by Kirchhoff (Ref. 12, p. 89), imply 
that on a free edge of a plate the shear force and the twisting moment are not neces- 
sarily zero. The physical significance of this apparent reduction in the necessary 
number of boundary conditions has been explained by Thomson and Tait (Ref.12, 
p. 90). 


3. Equilibrium Conditions 


Referring to Fig. 1, for moment equilibrium in a radial direction, 


In all the problems considered here the diameter of the hole is assumed large 
compared with the radial width of the reinforcement. This implies that the ¢/2 
term in expression (3) can be neglected, so that the equation becomes 


Similarly, for moment equilibrium in a direction tangential to the curve of 
the hole 


From Fig. 1 it is seen that the shear force acting on an elemental length of reinforce- 
ment does not act at the centre of its radial cross section. This shear force therefore 
gives rise to an offset moment at (Q,), 66/2, which tends to bend the reinforcement 
about a line tangential to the curve forming the hole. When dealing with particular 
problems, it can be verified that this moment is small compared with the tangential 
component of M and, therefore, can be neglected. 
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bent plate 


dS 

@Q.).=0. . ‘ ‘ 
It is worth noting that this relationship between the shear force and the applied load 
is the same as for a straight beam. 


From (4), (5) and (6), the following equations are obtained for the bending and 
twisting momenis acting on a radial section of the reinforcement : — 


@M 
“de +M=a 4 (M rata (M -@ (Q,)a. (7b) 

Writing (M,), etc. in terms of the deflection surface of the plate, the equations 
of equilibrium become 


de + =aD + * 
1 dw 
| 


(2-v) aw 3 dw dw 
oroe? or 


arLor + or * 


4. Compatibility Conditions 


For bending compatibility between the reinforcement and adjoining parts of the 
plate, the following boundary conditions must be satisfied :— 


(i) Ina plane perpendicular to the surface of the bent plate and which con- 
tains an elemental length of the reinforcement, the curvature of the 
reinforcement must be equal to the curvature of the plate. 


(ii) The rotation of a radial section of the reinforcement must be equal to the 
radial slope of the plate at the point considered. 


Since the radial width of the reinforcement is small compared with the diameter 
of the hole, an elemental length of the reinforcement can be considered as being 
straight when determining the moment necessary to produce its change in curvature. 


142 The Aeronautical Quarterly 


R. HICKS 
+ 
8a) 


PLATE WITH REINFORCED CIRCULAR 


Ficure 2. Curvature of reinforcement in a plane perpendicular to the surface of the plate. 


Thus, if M, is the component of the total moment M due to a change in curvature 


2 
M.= = El (ow 


where R, is the curvature of an elemental length of reinforcement in a plane 
perpendicular to the surface of the bent plate, as shown in Fig. 2. The negative 
sign has been inserted to be consistent with the convention that a positive moment 
bows the plate downwards. 


Referring to Fig. 3, at a depth z below its neutral axis the bending stress induced 
in the reinforcement due to its radial rotation is given by 


(10) 
where the negative sign has been inserted to be consistent with the convention 
previously stated. 


Using (10), the moment induced in the reinforcement due to its radial rotation 
is given by 


| 

a 

Et (2 a/2 
t{ow 
-d/2 
El (dw 
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Ficure 3. Section through deflected reinforcement and plate. 


Adding (9) and (12), the compatibility expression for the bending moment is 


low 1 

M=-El rot, oF (13) 

When deriving expression (13) it was assumed that the cross section of the 
reinforcement was rectangular; however, the relationship is true for a reinforcement 
having a cross section of any given shape, provided that the radial width of the 
reinforcement is small compared with the diameter of the hole. As for the bending 
moment, the twisting moment acting on a radial cross section of the reinforcement 
can be resolved into two components. The first component is proportional to the 
twist in the reinforcement and is independent of the diameter of the hole. The 
second component is proportional to the tangential slope of the reinforcement and is 
due to the fact that the twisted member is curved in a horizontal plane. 


Since the radial width of the reinforcement is small compared with the diameter 
of the hole, an elemental length of the reinforcement can be considered as being 
straight when determining the component of T due to the relative twist of two radial 
faces of the element. That is 


where K is a constant depending on the cross-sectional shape of the reinforcement 
and xr is the twist per unit length of reinforcement. 


( 
Now Xr=- 


EK (0?w 


The negative sign has been inserted because a positive twisting moment corresponds 
to a negative radial slope in the plate. 


Referring to Fig. 4, abcd is an infinitely thin elemental length of reinforcement 
lying in a horizontal plane. Assume that the elemental slice rotates about a 
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horizontal line oc so that the side ab is displaced by an amount dw to a’b’. Then 
the rotation of ab relative to cd is 'w/a. This implies that if an elemental length of 
reinforcement is given a tangential slope dw/(a5@) with no corresponding change in 
the radial slopes of its sides coniaining ab and cd, a twisting moment must be 
applied to the elemental length of reinforcement to suppress the rotation dw/a. 


Sw 


Ficure 4. Infinitely thin elemental length of reinforcement. 


Thus, for a unit length of reinforcement the suppressed twist is given by 


_ | ew 


That is, due to a tangential slope, the twisting moment induced in the 
reinforcement is 


EK ow 

Ts= EK\s= 06 . (16) 
From (15) and (16), the following expression is obtained for the twisting 

moment in the reinforcement :— 


1 1 ow 
dar 


T= -EK 


Values of K for reinforcements having cross sections of different shapes are given in 
standard works on strength of materials; however, it is worth noting that, for the 
important practical case of a reinforcement with a rectangular cross section, the 
appropriate expression for K is given by 


where d is the depth of the reinforcement and ¢ its radial width in the plane of 
the plate. 


It can be verified that the series part of expression (18) converges rapidly, so 
that K can be evaluated easily for any given value of t/d. 
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When the ratio d/t is sufficiently large, tanh (n~d/2t) can be taken as unity, so 
that the expression for K becomes 


Therefore, for a reinforcement with a square cross section, 


0:0703d* 
while an accurate evaluation of (18) gives 


0:06165d* 


Thus, in the extreme case, the approximate value of K represents an error of 
13-8 per cent. Since the accuracy of (19) rapidly increases as the ratio d/t increases, 
it follows that the approximate expression for K is sufficiently accurate for most 
rectangular reinforcements found in practice (say, when d > 21). 


5. Boundary Conditions 


Substituting (13) and (14) in (8), the boundary conditions around the hole 
reduce to 


1 1 1 1 dw (2- (2—v) 


6. Plate Deflection Surface 


The problem of a reinforced circular hole in a bent or twisted plate reduces, 
therefore, to finding a solution of 


7 
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which satisfies the boundary conditions (21) and the boundary conditions around 
the edges of the plate. 


When dealing with particular problems it is convenient to consider the 
deflection surface of the perforated plate as being composed of two parts as follows 


w=W,+W., ; . (24) 


where w, is the deflection surface of a similarly loaded unperforated plate and w, is 
the change in this deflection surface due to the introduction of the hole. 


For the particular problems considered here, the deflection surface of the 
unperforated plate is known, so that the solution of the problem is obtained by 
finding a complementary solution of (23) which satisfies the conditions around the 
edges of the plate and which, when added to w,, gives an expression for w which 
satisfies the boundary conditions around the hole. 


For a solution in polar co-ordinates, the deflection surface of the unperforated 
plate is given generally by the expression 


(1) cos S, (7) sin nb, 


where C, (r) and S, (r) are known functions of r. 


It is worth noting that various specific forms of (25) are given in Ref. 12, for 
example on pages 261, 268 and 293. 


When w, has the form given in expression (25), the complementary solution 
of (23), which can be made to satisfy the boundary conditions around the hole and 
give a close approximation to the boundary conditions around the edges of the plate, 
is given by 


A A, 


r r A, r n 
Wa = log + Ag? log” + (“2 + log”) cos 6+ cos + 


B, Bas na\.: 
+ (78 +B, rlog”)sino+ > (2 + sin (26) 


where A,;, B,, etc. are constants which can be determined by substituting 
w,+w, in (21) and solving the resulting equations by equating the coefficients of the 
cos n@ and sin n6 terms to zero. For convenience the arguments of the log terms 
have been made dimensionless by introducing the parameter s. When dealing with 
particular problems this parameter is given by a lateral dimension of the plate. As 
shown later for bent and twisted plates, differentiation of (26) results in expressions 
for the moment and shear components which do not contain the parameter s. 


If (26) is substituted in (1), expressions are obtained for the change in the 
moment and shear distribution in the plate due to the introduction of the hole. When 
dealing with practical examples, it can be verified that this change in the moment 
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and shear distribution is negligible except in regions of the plate immediately 
surrounding the hole. Thus, provided that the least lateral distance from the centre 
of the hole to the edge of the plate is several times the diameter of the hole, one need 
satisfy only the boundary conditions (21) when (26) is used to give the change in the 
deflection surface of the plate. 


This also follows from the work of Goodier''’ who considered the effect of the 
relative size of an unreinforced circular hole on the moment distribution in a bent 
plate. He found that when the ratio of the width of the plate to the diameter of the 
hole is reasonably large (say, greater than 6), the change in the moment distribution 
due to the introduction of the hole is local and independent of the relative size of 
the hole. 


7. Uniform Moment Distribution about One Axis 


Because the Principle of Superposition is applicable to plates having small 
deflections, the solution of this problem is particularly important, since it can be 
used to obtain solutions to various other problems of reinforced circular holes in 
bent or twisted plates. 


y 


Figure 5. Uniform bending x 
about one axis. 


Referring to Fig. 5, M, is the uniform bending moment per unit length of plate 
which produces bending about the xx axis. In this case it can be verified that the 
deflection surface of the unperforated plate is given by 


or, in polar co-ordinates, 


Comparing this with (26), and remembering that the boundary conditions (21) are 
satisfied by equating the coefficients of the cos n@ and sin n@ terms to zero, it follows 
that the change in the deflection surface is given by 

-M v A 23 


r 
| Aes lows + 20+ Ap, 608 28 | 
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It should be noted that there is no need to define the parameter s, since this 
does not appear in the resulting expressions for the moment and shear distribution 
in the plate and reinforcement. 


Substituting (27) and (28) in (1), it is found that 


M, 2 2 


6(1-¥) 


Ay, cos 26 — 


4 
A,,cos 26]. (29a) 


Me= v*)+2 (1 —v*) cos 26+ 


6(1-v) 


A,, cos 26 


2M. 


~ (i= 


When dealing with particular numerical examples, it can be verified that the 
terms containing the constants A,,, become relatively small and can be neglected 
when the ratio r/a is sufficiently large (say r/a > 3). This implies that the effect 
of the hole is local, so that the initial assumption concerning (26) is justified, that is, 
(26) can be used to determine the change in the moment and shear distribution if the 
constants A,,, afe determined by considering the boundary conditions around 
the hole. 


Referring to (29), it is seen that the shear forces Q, and Qs are due entirely to 
the presence of the hole. However, as previously explained, on the edges of the 
plate these shear forces are related to the twisting moment by the relationships (2), 
and these relationships must be used when considering the effect of the hole on the 
stresses at the edges of the plate. For example, for a similarly loaded annular 
plate, the shear on the boundary is given by 


M, 3(1-v) Az, (3-%) Ag, 


When r/a is sufficiently large, this is virtually the same as the shear on a similarly 
loaded unperforated plate. 


From (8), (13) and (17), the moments and shear force acting on a radial section 
of the reinforcement are given by the expressions (31) which follow. 
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M, El ig 
+2 (1+) 26+ ~ cos 29-4 cos 26] (31a) 
_M,_EK An 
T= Az) sin 26 
_ ~My Ax 


Adding (27) and (28), and substituting the resulting expression for w in (21), it 
is found that 


(32a) 
(1+) Ay —2(1+y) a? Ay=0. (32c) 


The solution of these equations gives 


Ay=-20-) 
= —(1+9) (3+ 7 (1+%)—A (3—v)—yA (1 (33) 


+A G+ —9)] 


(3-(1 —vA+y)—yA (1 a? ie 


For a similarly loaded plate with an unreinforced circular hole y=A=00 and 
the expressions (33) reduce to 


A,;=2 a? 
Ay,» (34) 


Substituting these values for the constants in (29), it is found that the maximum 
moment in the plate acts in a tangential direction and occurs on the boundary of 
the hole at the point 6=0. Furthermore, this maximum bending moment is 
given by 


§+3y 
— 


Thus, when v=0-3, the effect of the hole is to increase the maximum stress in 
the plate by 74:2 per cent. 
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APPLIED MOMENT M 


APPLIED MOMENT My 


Ficure 6. Numerical example showing variation of the bending moments along the radial 
lineg 6=0 and 6=7/2. 


8. Numerical Example 


To illustrate the effect of a reinforced hole on the moment distribution in a 
plate, consider the numerical example where the plate and reinforcement have the 
following dimensions : — 


Thickness of plate =h 
Depth of reinforcement =3-56h 
Radial width of reinforcement =h 
Diameter of hole =41h 
Poisson’s ratio =0-3. 


For this particular example, Fig. 6 shows the variation in the radial and 
tangential bending moments along the radial lines 6=0 and 6=/2. As would be 
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ry 


Ficure 7. Radial and tangential bending moments in the plate around the hole at the 
point 6=7/2. 


Fioure 8. Radial bending moment in the plate around the hole at the point 6=0. 
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Ficure 9. Tangential bending moment in the plate around the hole at the voint 6=0. 


expected, the change in the stress distribution is confined to a small area immedi- 
ately surrounding the hole. That is, the solution for an infinite plate is applicable 
to a plate of finite size, provided that the shortest distance from the centre of the 
hole to the edge of the plate is several times greater than the diameter of the hole. 


This example illustrates that a suitable reinforcement results in a complete 
reduction in the stress concentration, so that the stress in the plate with the 
reinforced hole is nowhere greater than in a similarly loaded unperforated plate. 


9. Design Data 


For design purposes the choice of a suitable reinforcement is simplified if the 
moment concentrations at discrete points around the hole are known for reinforce- 
ments having various cross sections. For this reason, values for radial and 
tangential moment concentrations have been evaluated at the points #=0 and 6=7/2 
for a wide range of reinforcements with different bending and twisting stiffnesses. 
These values are based on the assumption that Poisson’s ratio is 0:3, and are given 
by a series of curves shown in Figs. 7, 8 and 9. 


Since the moment concentration (i.e. stress concentration) around an unrein- 
forced hole is 1-74, it follows from these curves that reinforcements with excessive 
bending and twisting stiffnesses can produce stress concentrations greater than those 
in a similarly loaded plate with an unreinforced hole. 


It is worth noting that, because the Principle of Superposition is applicable to 
bent and twisted plates, the foregoing design data can also be used to determine the 
corresponding stress concentrations in plates under the action of various other types 
of applied moment distributions. 


10. Uniform Moment Distribution about Two Axes 


If M, and M, are the uniform moments per unit length of the plate which 
cause bending about the yy and xx axes respectively, the solution to the previous 
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problem is still valid provided that the following changes are made in the formulae 
for the deflection surface and bending moments, and so on. 


(i) For terms which do not contain 4, substitute (M,+M.) for M,. 
(ii) For terms containing @, substitute (M, —M.) for M,. 
Thus, as an example, (29a) would become 
(M,+Ms,) (M,—M:) 


6(1 —v) Ags + 


} cos 20. 
The expressions for the constants A,,, A,;, and A,, remain unchanged. 


11. Uniform Twisting about Two Axes 


If M; is the uniform applied moment which twists the plate about the x and y 
axes, it can be verified that 


Myxy M,r’ sin 26 


2(i-v)D ’ (37) 


+ B,,) sin 28, 
(3+ (1 —v)"] 
[3+ (A+ ¥) + yA (3 + — 
[3+ (A+ 75+ ¥)+ yA (3+ —v)] 


B,, 


My 3 
20-9 Bu sin 26 


{1-72 cos 20 


Ba 


By 
Qe= - —* cos 20 


so that 
where 
Also 
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D(i-») 


M,EK 3B,, ” 


B,, 


S= 


12. Non-Uniform Applied Bending and Twisting about x and y Axes 


When the applied bending and twisting is non-uniform about the x and y 
co-ordinate axes, the deflection surface of a similarly loaded unperforated plate is 
given by 


cos6+ B,.r° sin6+ A,.r"*?) cos nO + 
n=2 


+ = (Bur + B,,r***) sin . . (41) 


It can be verified that the A,,r° cos @ and B,,r° sin @ terms individually satisfy 
the boundary conditions (21), so that the complementary constants, A,;, A,,, Bis 
and B,,, in (26) are all zero. Furthermore, since the boundary conditions (21) 
involve an even number of differentiations of w with respect to @, it follows that for 
any value of n the relationships between the constants of the cosine terms (that is, 
Any» Anz, Ans and A,,) ate the same as the relationships between the constants of the 
sine terms (that is, B,,, B,., Ba; and B,,). 


Adding (26) and (41), and substituting the resulting expression for w in (21), 
around the hole the boundary conditions are satisfied when 


2Ao, [1—y (1+ Aos (1 + (1 —»)Ja-? =0 (42a) 


Ayn [(n—1)—A (1 a*-? + Ans {n (n+ 1)—A [n (1 + 4] } 
— [(n+ 1)+A (1 —A,, (42d) 
[((n—-1) -y (1 + Ang 1)—-y [n (1 —v) (1+) ]} + 
+ Aggn [(n+ 1) +7 (1 + An, + + 
(42c) 


where the equations defining the B,,, constants are the same as these equations 
defining the A,», constants. 
Solving the equations (42), it is found that 
[1 (1+ A,,a’ 
{1+y(1—v)] 
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—2 


A,sna-**{(n? — 1) + (A+ yX2n+ 1+ v)+ yA (3+ v)} 
=A,,n(n—1){(n? —1)+y (1+) —A (3 —v) —yA (1 —v)?} + {n? (n? — 1) + 
+ ny (3+ ¥)—A[(S—v) n? —8]—yA [n? +8 (1+ )]} . (435) 


Ana~?*{(n? — 1) + (A+ y2n+ 1+9)+7A (3+ —v)} 
= —A,,na~?{(n? —1)—(1 — vy + A) — — ¥)?} — Ans (n+ I) {(n? — 1) + 
+y(1+v)—A (3 (1 v)*}. (43¢) 


When the deflection surface for a bent or twisted unperforated plate is known, 
the expressions (43) completely determine the unknown constants in the general 
expression for the deflection surface of a similarly bent or twisted plate with a 
reinforced circular hole, it being understood that identical expressions define the 
Bam Constants. 


Substitution from expressions (43) into equations (1) and (2) gives the following 
general expressions for the moments and shear forces in the plate and 
reinforcement : — 


{2(1+¥) Ay, —2D (3 + v) A cos 6 — 


-D = {n(n—1\1 + [n? —v) +2 (3+ ¥)+2 (1+)] + 


+n(n+1\(1—v) Ot? + 
+[n? cosn@ . (44a) 
Me= —D {2(1+%) Ay, + (1 A, — 2D (1 + 39) A,r cos 6+ 
+D = {n(n—1)1 —») + [n? (Ll —v) (1 + 39) —2 (1+ 9)) + 
+n(n+1X1—v) + 
+ [n? (1—v)+n(1+3v)—2 (1+ %)] Aar-"} cosné6 . (44d) 
—2(1—¥) DA, r n[(a—1) + (n+ 1) Aur — 
—(n+1) +” —(n—- 1) sin nO 
—4DA,,r~-'—8DA,, cos 6— 
-4D = n[(n+ 1) 1) cos . (44d) 


Q,=8DA,, siné6+4D [(n+ 1) An —(n— 1) sin (44e) 
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= — EI {2A,, + —2El A,,a cos 0+ 
+El & [n(n—-1) (n+ + 
n=2 


+n(n+1) A,a-“+” + (n—1)(n+ 2) cos 


T=EK A,,asin9+ EK [(n— 1) 


+(n+ 1) A,,a" —(n+ 1) —(n— 1) A,,a~"] sin n6 


S=8DA,,.a sin 1) A,,a" +(n—1) A,,a~-"] sin nO + A, 


where A is a constant of integration. 


It is outside the scope of this paper to consider in detail the application of the 
solution for non-uniform applied bending and twisting about the x and y co-ordinate 
axes. However, as an example, reference can be made to the case of a reinforced 
circular hole in a simply-supported equilateral triangular plate bent by constant 
moments distributed around its boundary. In this case, the deflection surface of a 
similarly loaded unperforated plate is the same as that of a uniformly stretched and 
uniformly loaded membrane, and is given by expression (45), which was first obtained 
by S. Woinoasky-Knieger (Ref. 12, p. 102). 


M, 


where M, is the moment per unit length applied to the boundary of the plate. 


Therefore, for this problem, 


and all other constants in equation (41) are zero. 


Thus from equations (43), (44) and (46), expressions can be obtained for the 
moments and shear distributions in a plate with a reinforced hole. 
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Experimental Investigation of the Stalled Flow 
of a Single-Stage Axial Flow Compressor 


M. D. WOOD, J. H. HORLOCK* and E. K. ARMSTRONG+ 
(University Engineering Laboratory, Cambridge) 


SumMaRY: Experimentally it is shown that variation of axial spacing between 

neighbouring blade rows can alter the stall characteristics of a compressor 

stage. The stall cell patterns of this compressor stage are described in detail 

and changes of slope in stage temperature and pressure rise characteristics are 

explained in terms of stall propagation. Strain gauge tests on stator blades 

are used to indicate the possibility of large dynamic stresses occurring when 
the stall cell frequency is equal to the blade natural bending frequency. 


1. Introduction 


A paper by Horlock‘” has described theoretical and experimental investigations 
of the “interference” effect of one blade row upon the performance of neighbouring 
rows. The investigations concerned time average effects in an axially symmetric 
field of flow. In extending this work from the normal to the stalled flow region it 
was decided to attempt to gain some insight into the mechanism by which each blade 
row Stalled. Other experimenters” had reported non-uniform flow during the first 
phase of stall, and this phenomenon had become known as “Stall Cell Propagation.” 
It was found that sharply defined areas of high loss rotated in the blade rows at 
angular velocities ranging from 25 to 65 per cent of rotor velocity. Accordingly, it 
was decided to use a two-channel hot-wire apparatus developed in the University 
Engineering Laboratory to determine whether the axial flow compressor described 
in Horlock’s paper was subject to stall propagation. After confirmation of the 
qualitative results of Ref. 2, the effect of blade row interference on the appearance 
and frequency of the stall cell patterns of a single stage compressor was investigated. 
This paper lists the number and frequency of stall cells at different flow points for 
an isolated rotor row, an isolated stator row, a rotor-stator combination, and finally, 
the complete stage consisting of guide-vane, rotor, and stator rows. Since the main 
interest of the compressor designer in stall propagation concerns the possibility of 
the stall cell exciting blades at a resonant frequency, strain gauges were fitted to a 
stator blade and a measure of the dynamic stresses occurring under these conditions 
was recorded. The last section of this paper compares the stress levels which were 
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measured in the normal aluminium stator blades with those in specially designed 
plastic blades, whose natural frequency in bending was close to stall cell frequency. 


This work was done at Cambridge University Engineering Laboratory during 
the years 1953-1955. 


NOTATION 
C., mean axial velocity 
U,, blade speed at design radius 
K, specific heat of air at constant pressure 
AT stage total temperature rise 
AP stage total pressure rise 
air density. 


> 


2. Apparatus and Instrumentation 
2.1. COMPRESSOR 


The single-stage compressor was driven at 6,000 r.p.m. through a 2-086 :1 step- 
up gear box by a 25 h.p. induction motor. The rotor tip speed was 370 ft./sec. 
Provision was made for varying the axial spacing of the blade rows. The complete 
stage of guide vane-rotor-stator had been designed to simulate the first stage of a 
conventional aircraft compressor. In the design the exit axial velocity and rotor 
absolute air outlet angle were assumed to be constant at all radii. The mass flow was 
controlled by an outlet throttle and measured by reference to the static pressure 
depression in the inlet ducting. The temperature rise was deduced from the torque 
input, allowance being made for friction and windage. A mean delivery total 
pressure was obtained by adding an exit dynamic head to the mean of the inner and 
outer wall static pressures. The mean exit dynamic head was computed on the 
assumption that the flow through the compressor was incompressible and that the 
angle of flow in the outlet plane was equal to the design value. Fig. 1 shows 
the general arrangement of the compressor. 


Details of the blading are as follows :— 


blade tip radius 7-0 in. 
blade tip speed 370 ft./sec. 
blade hub radius 2°8 in. 
hub/tip ratio 0-4 

number of rotor blades 19 

number of stator blades 23 

number of inlet guide vanes 23 

blade chord 1:25-1°5 in. 


C5 profiles on circular arc camber lines. 
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Blade angles: 
Rotor tip inlet 57°, outlet 42°, space/chord ratio 1-50, 
hub inlet 25°, outlet 4°, space/chord ratio 0-60. 
Stator tip inlet 44°, outlet 24°, space/chord ratio 1-28, 
hub inlet 28-5°, outlet 3-5°, space/chord ratio 0-60. 
Inlet guide vanes :— 
tip outlet 44°, space/chord ratio 1-39, 
hub outlet 10°5°, space/chord ratio 0-60. 
Axial distance from the trailing edge of a blade row to the leading edge of the 
adjacent downstream blade row (rows closely packed) is 0-69 in. 
Axial distance from trailing edge of second stator, when fitted, to leading edge 
of outlet guide vanes is 5-0 in. 


2.2. Hot-WiRE MEASUREMENTS 

The fluctuating velocities and stall cells were investigated by means of two 
constant-current hot-wire anemometers, feeding, through amplifiers, the two inputs 
of a double beam oscilloscope. (The frequency characteristics of the amplifiers com- 
pensated for the reduced response of the hot wires to high-frequency velocity 
fluctuations.) When stall cell propagation is established in a compressor, the 
frequency of excitation of a fixed blade may be determined with a single anemo- 
meter. To determine the velocity of propagation of a single cell, and hence the 
number of cells, the time lag between signals from two anemometers at several 
known angular positions must be measured. The double beam oscilloscope offers 


Ficure 1. General arrangement of axial flow compressor. 
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Ficure 2. Hot-wire anemometers arranged for measurement of stall propagation velocity in 
compressor blade row. 


a convenient method of performing these measurements. In practice a single sweep 
of the two oscilloscope traces was photographed at each of four angular spacings of 
the nrobes covering, in 9° steps, the range from 99° to 126°. Radial traverses of a 
single proble showed the extent to which the stall cells covered the annulus. Fig. 2 
shows the compressor with two probes in position for a stall cell investigation. 


2.3. BLADE STRESSES 


Strain gauges were mounted on the concave surface of the stator blades at right 
angles to the chord at the aerodynamic tip.* Since the blades were cantilevered 
from the outer casing of the compressor, the gauges measured the maximum strain 
due to bending in the fundamental mode. The amplified out-of-balance signal from 
a Wheatstone bridge, of which the strain gauge formed one arm, was displayed on 
one cathode ray tube of a multi-tube recording unit. A hot-wire trace and a timing 
signal were shown on other tubes and the three were photographed simultaneously 
on a moving film. The strain gauge circuit was calibrated on each trace, thus 
allowing the blade stress to be computed from the known gauge factor and from 
Young’s modulus. The value of Young’s modulus for the plastic material was 
determined by comparative tip deflection tests between the geometrically identical 
plastic and aluminium blades. The value of Young’s modulus for the aluminium 
was known. 


*Throughout this paper, root and tip refer to the aerodynamic root, inner radius, and the 
aerodynamic tip, outer radius, respectively. 
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Ficure 3. Stall propagation on two-dimensional cascade. 


3. Qualitative Explanation of Stall Cell Propagation 


A simple explanation of stall propagation has been advanced by several writers. 
It is included here for completeness. 


If a cascade, as shown in Fig. 3, is at an incidence close to, but below, the 
stalling value, one blade (A) may, due to some misalignment, fault in manufacture, 
or damage due to a number of possible causes, stall before the rest. This results in 
a blockage of the passage AB and some of the inlet flow will be diverted to neigh- 
bouring channels BD and AC. Clearly, the effective incidence of blade B will be 
increased and that of C decreased. Blade B may then stall, and the disturbance 
therefore propagates in the direction A to B to D. As the region of stall grows, 
blade A is at a position where, by the previous argument, the incidence is now 
decreased. At some time after the initial triggering action, the incidence on A is 
sufficiently reduced for it to unstall again. The disturbance then continues to 
propagate as a “cell” with well defined boundaries. 


4. Compressor Arrangement and Stall Patterns 


The compressor may be fitted with a maximum of five closely spaced rows of 
blades (followed in all cases by a row of guide vanes downstream in the outlet duct). 
The five rows comprise guide vanes, first rotor, first stator, second rotor and second 
stator. This arrangement would be designated GRSRS/G. Subject to the limita- 
tion of there being at least one rotor, any or all of the other rows may be omitted. 
Any axial interference effect between blade rows separated by at least two blank 
rows should be small and the elements of a single stage may therefore be assembled 
as follows :— 


Isolated rotor Represented by G- -R-—/G 
Isolated stator Represented by GR- —S/G 
(Guide vane-rotor combination also represented by GR-— —S/G) 
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Rotor-stator combination Represented by G— —-RS/G 
Simple stage Represented by GRS- —/G. 


The stall patterns associated with the various flow points of the characteristic 
of each compressor element are shown diagrammatically in Figs. 4-7. To simplify 
description of the patterns a notation has been developed which is best explained 
by an example. In Fig. 4, at a flow represented by C,,/Un=0°492, the pattern is 
described by RT/2/67/134. The sketch shows that stall propagation is occurring 
at the rotor tip (RT) and that 2 cells are rotating in the annulus. The designation 
RT /2/67/134 supplies the extra information that each cell is rotating at 67 revolu- 
tions per second and that therefore the fundamental excitation frequency, relative 
to a blade fixed in space, is 134 c./sec. Similarly, SH (E)/ 1/86/86 indicates that one 
cell is propagating on the stator hub in an erratic manner at approximately 86 r.p.s. 


In noting the erratic propagation at this flow it must be emphasised that all the 
stall cells observed showed some degree of unsteadiness. Both the frequency of 
appearance and extent of the cells varied in a non-systematic manner with time. The 
values given are mean values obtained by averaging the frequency and extent of 
several cells selected by photographing oscilloscope traces at random intervals of 
time. No attempt was made to measure the degree of unsteadiness of the stall cells 
nor the magnitude of the velocity fluctuations forming the stall cells. 


HOT WIRE TRAVERSE 


---4040 


4 0-30 


0-20 


0-10 


KEY 
@Q—STALL CELLS 
—-ERRATIC STALL CELL 
PROPAGATION 
STALL CELLS FROM 
PREVIOUS ROW 


E -—-ERRATIC PROPAGATION 


——CONTINOUS STALL 
AROUND THE ANNULUS 


IRT/3, 66/198 


Ficure 4. Single-stage compressor. Stall cell patterns. Compressor build G - - R - /G. 
164 The Aeronautical Quarterly 


STALLED FLOW IN COMPRESSORS 


5. Experimental Results 
5.1. STALL PATTERNS 
(a) G--R-/G 

Figure 4 shows the results obtained for the isolated rotor. The method of 
evaluating the outlet total pressure outlined in Section 2.1 is of limited value when, 
as in this compressor arrangement, there are large whirl velocities in the outlet plane. 
For this reason the total pressure rise characteristic is not shown in Fig. 4. 


All the characteristic phenomena associated with stall propagation are 
demonstrated in this simple compressor. Apart from slight thickening of the blade 
wakes shed from the rotor, no disturbance was noted as the flow was reduced until 
the point C,,/U=0-503 was reached. At this point a marked change in slope 
occurred in the characteristic. The hot wire showed the sudden appearance of a 
single propagating stall cell on the rotor tip. Further reduction in flow caused the 
stall cell to split into two, three and later four cells. This is characteristic in that 
a reduction in flow tends to increase the number of cells rather than to increase the 
size of the original single cell. Setting the flow to a value which lay between two 
regions of stable propagation resulted (as at C,,/U,=0-438) in the flow oscillating 
in an arbitrary manner between the patterns corresponding to slightly greater and 
slightly smaller flows (three- and four-cell configurations respectively). 


The change in slope of K,AT/U,,? at C,,/Um=0°503 is explained by a decrease 
in the fluid deflection associated with the appearance of a stall cell. The second 
change in slope, at C,, /U,=0°495, is not yet clearly understood. It is thought that 
the increase in temperature rise at low flows may be a measure of the energy which 
is absorbed in recirculation of the flow at the blade tips. 


At flows below C,,/U,=0-359 the three-cell pattern degenerated into a band 
of continuous stall in which all rotor blade tips were stalled at all times. 


(6) GR--S/G (Fig. 5) 

The isolated stator showed propagation at a higher flow than on the isolated 
rotor. Two cells were found which propagated in an erratic manner on the stator 
hub. The presence of these cells is associated with a distinct flattening of the 
pressure rise characteristic. 


Stall propagation was observed on the stator over a wide range of flows, the 
interference from the rotor tip stall being limited to a reduction in the amplitude of 
the stator stall cell disturbances. “Shadows” of rotor tip stall cells were found 
downstream of the stator tip. They appeared as small flow disturbances 
synchronised with the rotor stall patterns. 


This same compressor arrangement demonstrates the effect of placing a guide 
vane row immediately upstream of the rotor. The single- and two-cell configura- 
tions appeared at the rotor tip at substantially the same flow as before, but the 
three- and four-cell patterns were not observed. The two-cell patterns persisted as 
the flow was reduced until, at C,,/U,,=0°443, the rotor tip stalled continuously. 
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The rotor hub was unstalled and at mid-blade height erratic propagation was 
observed. Further flow reduction resulted in an increased area of complete tip stall 
with consequent reduction in the hub area of normal flow. As before, the inter- 
mediate region showed erratic propagation. 


(c) G——RS/G (Fig. 6) 


The rotor-stator combination of this compressor arrangement showed stall 
patterns which were markedly different from those for each of the rows in isolation. 
As the flow was reduced the stator hub stalled at C,, /Um=0-515 (approximately the 
same value as for the isolated stator). As a result of this stall, the rotor tip flow was 
relieved to such an extent that stall, defined by the discontinuity in slope of the 
temperature rise characteristic, was delayed until C,,/Un=0°455. At this flow a 
three-cell pattern was formed on the rotor tip, thus in turn relieving the stator hub 
and causing it to unstall. Further reduction in flow produced continuous rotor tip 
stall and, at a very low flow, erratic propagation was observed on the stator hub 
once more. 


As the throttle was opened again marked “hysteresis” was observed. The rotor 
continuous stall broke down into a single propagating cell which persisted until an 
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increasing flow of C,,/U,=0°440 was reached. This explains the loops in both 
the temperature rise and pressure rise characteristics. It must be emphasised that 
this “hysteresis of pattern,” in which different stall cell configurations are observed 
in decreasing and increasing flow, is inherently different from the “hysteresis of 
appearance”, in which the same configurations appear at slightly different flows, 
reported by other experimenters. With the single exception of this compressor 
arrangement (G— —RS/G), no other hysteresis was observed. In all other cases 
stall patterns were established, within the limits of experimental error, at the same 
values of C,,/U,,, whether the flow was increasing or decreasing. 


(d) GRS- —/G (Fig. 7) 


This arrangement represented the complete compressor stage. As the flow was 
reduced, stall propagation was observed at the stator hub at approximately the same 
flow as in the isolated stator (C,,/U,=0:530). The four-cell pattern, which was 
observed at C,,/U,=0°530, degenerated into a very erratic two-cell pattern as the 
flow was reduced further. The interference effect of the stator hub stall delayed 
the rotor tip stall until C,,/U,=0-450 was reached. At this flow one large rotor 
tip cell was formed, which, rather surprisingly, did not unstall the stator. However, 
the flow at the stator root was so disturbed that it was not possible to determine the 
number of cells. Further throttling caused the rotor tip cell to split into two cells 
(C., /Um=0-408) before joining up to form a continuous band of stall. When this 
band was established (C,., /U,, =0°309) the interference effect was sufficient to unstall 
the stator hub, which only restalled at very low flows (C,, /Um=0-202). 
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As in the arrangement G— —RS/G, shadows from the rotor tip stall were 
observed behind the stator tip. 
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5.2. DYNAMIC STRESSES IN STATOR BLADES 


It was felt that compressor blades which were subject to periodic flow changes 
at their resonant frequency might show large resultant stresses. In order to investi- 
gate this possibility, a strain gauge was mounted on one of the aluminium stator 
blades of the compressor arrangement GRS— —/G. On the opposite side of the 
annulus, six of the aluminium stator blades (first natural frequency in bending 
375 c./sec.) were replaced by plastic blades (natural frequency 120 c./sec.) of 
identical geometric form. One of these plastic blades was also fitted with a strain 
gauge. It was found that the replacement of aluminium by plastic blades in no way 
altered the aerodynamic performance of the stage. 


In the normal running condition the mean dynamic stress in the aluminium 
blade was much higher than that in the plastic blade. The stress reduction in the 
plastic blade was considered to be the result of the high material damping. The 
logarithmic decrement of a free vibration of the plastic blade in its mounting was 
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found to be approximately one hundred times greater than that for the aluminium 
blade under the same conditions. To form a basis of comparison, all dynamic blade 
stresses were expressed as multiples of the mean peak-to-peak dynamic stress of the 
same blade under design flow conditions. The “Amplification Ratio” is shown for 
the plastic and aluminium blades at several compressor flow points (Fig. 8). 


At C,,/Um=0°513 the stator hub is subject to erratic stall propagation of 
approximate frequency 80 c./sec. The stress amplification for the aluminium blade 
(natural frequency 375 c./sec.) is less than 2, while that for the plastic blade (natural 
frequency 120 c./sec.) is 10. 


At C,,/Um=0°462 the mean value of the stall propagation frequency had 
increased to approximately 100 c./sec. The magnitude of the velocity fluctuations 
forming each stall cell was slightly greater. This resulted in the amplification ratio 
of the aluminium blade increasing to 4, while that of the plastic blade reached its 
maximum value of 45. The aluminium-blade strain-gauge record showed a small 
but clearly defined increase of dynamic strain at approximately four-cycle intervals. 
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Evidently excitation of the aluminium blade by the fourth harmonic of the stall cell 
frequency was not serious. However, the increase of stall cell frequency to a value 


more nearly equal to the natural frequency of the plastic blade resulted in large 
siress amplification. 


At C,,/Um=0:309 the rotor tip was completely stalled and the stator hub 
completely unstalled. The peaks which had been observed on the strain-gauge 
record of the aluminium blade at approximately four-cycle intervals at C,,/Un= 
0:462 were no longer observed. The aluminium blade amplification ratio had 
reached a value of 5, while that corresponding to the plastic blade had decreased to 
18. This decrease is attributed to the absence of excitation at a frequency near 
resonance. 


In these experiments the absolute magnitudes of the dynamic stresses were not 
excessive. Only the largest instantaneous values in the aluminium blade were equal 
to the mean aerodynamic stress, and, due to higher damping, the plastic blade 
stresses were much lower than those in the aluminium. 


6. Conclusions 


Marked interference effects were demonstrated on the performance of blade 
rows which were closely spaced, as in the conventional multi-stage axial compressor. 
The stall of one row may delay the stall of a neighbouring row or even, in extreme 
cases, unstall a previously stalled row. 


In all cases investigated, the flattening of the compressor pressure-rise 
characteristic was associated with stall propagation in the stator blade row, and 
discontinuities in the temperature rise characteristic with stall propagation in the 
rotor row. 


Amplification of dynamic stresses in compressor blades was clearly demon- 
strated when the fundamental stall cell frequency approached the natural bending 
frequency of the blade. In the single case in this investigation where a blade’s 
natural bending frequency was a multiple of the stall cell frequency, the strain- 
gauge record showed evidence of the effect of the periodic disturbance, but the 
resultant amplification of dynamic stress was negligible. 


In the various blade row arrangements of this one simple constant-speed 
compressor, a wide range of fundamental stall cell frequencies was encountered. A 
stator blade would have been subjected to periodic disturbances at 44, 62, 70, 85, 
100, 124, 135, 195 and 270 c./sec. It is therefore concluded that there is an urgent 
need for a precise theoretical treatment of the problem, in order that the appearance 
and range of excitation frequencies of the stall cells may be predicted from 
compressor design data. 


REFERENCES 


1. Hortocx, J. H. Experimental and Theoretical Investigations of the Flow of Air 
Through Two Single Stage Compressors. R. & M. 3031, 1955. 


2. Emmons, H. W., Pearson, C. E. and Grant, H. P. Compressor Surge and Stall Propaga- 
tion. Paper No. 53-A-65 A.S.M.E. Annual Meeting, November 1955. 


170 The Aeronautical Quarterly 
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Roughness from Insects 
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SUMMARY: In Reference 3, attention is drawn to the difficulties of measuring 
the streamwise extent of the roughness from insects. The present paper deals 
with the problem theoretically for an aerofoil in two-dimensional, incompres- 
sible flow. A tentative approach to the determination of effective excrescence 
height downstream of the leading-edge zone is also advanced. The application 
of these investigations, in conjunction with the analysis from Ref. 3 regarding 
the critical conditions for premature transition, leads to estimates of the 
amount of significant roughness which are in good agreement with flight 
observation. 


1. Introduction 


The roughness produced by insects during flight at low altitudes has become 
a familiar and significant problem with the development of the low-drag aerofoil and 
the full-scale investigation of boundary-layer control. Some preliminary experiments 
to measure the streamwise extent of the resultant deposits under various flying 
conditions are reported in Refs. 1 and 2. Certain wind-tunnel experiments® have 
also been performed with live insects to examine in more detail the geometrical 
properties of this natural form of roughness and to establish the conditions under 
which the flow will be made turbulent if steps are not taken to prevent contamination 
in the region where it is significant. At the same time, as explained in Section 3 of 
Ref. 3, the relevant technique of injecting insects into the flow at a point upstream 
of the model does not reproduce free-air conditions sufficiently well to give a strictly 
representative measure of the chordwise extent of the roughness. On the other hand, 
flight experiments designed to cover adequately, and independently, the principal 
variables involved are by no means devoid of complexity and, in addition, are 
inherently costly and laborious to perform. However, the streamwise extent of the 
roughness may be determined theoretically and, although the solution necessitates 
step-by-step calculations for specific initial conditions, it has seemed worth while, in 
view of the experimental difficulties just mentioned, to investigate the case of two- 
dimensional, incompressible flow past an aerofoil of arbitrary shape. This work is 
summarised in the present paper. 
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A simple theory is also advanced to establish the envelope of excrescence height 
(as defined in Section 4 of Ref. 3) within the predominant region containing deposits 
of an essentially fluid character (haemolymph), i.e., excluding the heavily roughened 
section immediately adjacent to the leading edge where the excrescences consist, 
more or less, of whole insects. Thus, since the chordwise extent of the roughness 
and the envelope of excrescence height are known, the extent of the significant 
roughness follows immediately from a knowledge of the critical conditions for 
transition. The relevant criterion has been established in Ref. 3. 


NOTATION 


O’x’,O’v’ co-ordinates of the particle P, taken at rest with respect to the 
undisturbed stream (Fig. 1) 


x’,,), initial position of P relative to the aerofoil origin O 


Ox,,Oy, aerofoil co-ordinates, with Ox, in the direction of motion of the 
aerofoil (Fig. 1) 


Ox,Oy aerofoil co-ordinates, with Ox measured along the chord of the 
aerofoil (Fig. 1) 


t time 
U, velocity of the aerofoil relative to O’x’, O’y’ 
W induced velocity of the air relative to O’x’, O’y’ (Fig. 1) 
U component of W in the direction of O’x’ 
V_ component of W in the direction of O’y’ 
y=tan-'(V/U) 
Q=W/U, 
w velocity of P relative to O’x’, O’y’ (Fig. 1) 
@ angle which the tangent to the particle path makes with O’x’ (Fig. 1) 
o=w/U, 
q={W? —2Ww cos 9) + w?}"/? 
c aerofoil chord 
a angle of incidence 
A=U,t/c 
| body length of insect 
4 volume of the insect 
p, density of the insect 


p air density 
D insect drag force 


Cy drag coefficient of the insect based on its relative air velocity and 
the area 


k=Cypp/(2lp,)_ (Fig. 3) 
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A,,A,,0 defined by the expressions in (7) 
x;y, cO-ordinates, with respect to Ox, Oy, of a point of impact 
X,, yr Values of x,, y, at the limit of roughness 
X;, chordwise extent of the roughness, measured from the leading edge 
5x, = 
X,; Value of x, on the Joukowski aerofoil 
Wy resultant velocity of P (Fig. 2) relative to Ox,, Oy,, with w, =w,/U, 
o angle which the tangent to the particle path makes with Ox, (Fig. 2) 
(i), Og (tr) values of wy at x;, and x,, y, respectively 
o(i),o(r) values of o at x, y, and x,, y, respectively 


w, impact velocity (component of wg normal to the surface), with 
= W,/ U, 


w, rupture velocity of an insect, with o,=w,/U, 


(w;) min Minimum rupture velocity of the insect aerial population with 
min (Wy) min / U, 


t slope of the profile with respect to Ox, Oy 
z(r) value of 7 at x,, y, 
7 (X13) \ values of 7 and d7/dx for an aerofoil of any shape 
(dz /dx);,, | at the section x= x,, 
& effective excrescence height 
n co-ordinate normal to the aerofoil contour 


w, velocity of the mass centre of the particle normal to the aerofoil 
contour 


n, distance of the mass centre of the particle from the surface at the 
instant of particle contact (w, =w,) 


n, distance of the mass centre of the particle from the surface ‘after 
completion of impact (w, =0) 


d=n,-n, 
d, value of d when w, =, initially 


Suffix 
J denotes quantities which relate to the Joukowski aerofoil 


2. Insect Motion in the Induced Velocity Field 


Before an analysis can be made of the conditions at impact, it is necessary to 
consider how the free activity, or random motion, of the insect in flight may affect 
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the problem. This factor is discussed fully in Appendix I, where it is shown, from 
the energy and drag data of Ref. 5, that, for the small, weak-flying species relevant 
to the present investigation’, the propulsive force which can be exerted is negligible 
compared with the drag force generated by the induced velocity field adjacent to the 
aerofoil. Consequently, in the problem under discussion, the significant insects may 
be regarded dynamically as inert particles, and will be treated as such in the present 
section, apart from the assumption that they continuously generate lift sufficient to 
react their weight. 


Ficure 1. 


Take, then, rectangular co-ordinate axes O’x’, O’y’ (Fig. 1) which are at rest 
relative to the undisturbed air. At time t=0, the origin O of the aerofoil is con- 
sidered to be at some large distance x’, from O’ such that the induced velocity due 
to the approach of the aerofoil with uniform speed U, along O’x’ is then vanishingly 
small*. Simultaneously, let a particle (insect) P be located on O’y’ at y’,. Hence, 
for any subsequent time f, the aerofoil origin O will be displaced x’,+U, t from O’, 
and the new co-ordinates x’, y’ of P will be 


t t 
x= | woos dt, yay y+ | wsing de, (1) 
e 


where the particle velocity vector w is defined (see Fig. 1) by w=we'*. Correspon- 
dingly, let W=We be the induced velocity produced by the aerofoil, the 
calculation of this quantity being considered in Appendix II. For the present, 
therefore, we assume that W is known, and proceed to obtain the hodograph of the 
particle. As shown in Ref. 5, the force exerted by the fluid on the particle may be 
approximated to the drag D, which will be written as 


(w- wy? { Ww cos w+iW sin ol, (2) 


*This condition is discussed in Section 3, where it is shown that x’, = —10c (c being the wing 
chord) is a satisfactory initial value. 
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where Cp, /, are, respectively, the drag coefficient, body length and volume of the 
insect, p is the air density and 


q= {W? —2Ww cos — 9) + w?}?/?, 


Hence, separating real and imaginary parts of (2), we obtain immediately for the 
equations of the hodograph 


dw 
= =kq{W cos 


wit =kqW sin 


in which k is a coefficient having the value 


and m,, p, are the mass and density of the insect respectively. It will be convenient, 
however, in the subsequent analysis to write 


Q=W/U,, w=w/U,, A=U,t/c, «=Cypc/(2pd. 
In this case, equations (3) and (4) may be expressed non-dimensionally as 


=KA, (Q-w), 


9). 


® 


where 


2PQe 
{1+ 


ra 


Approximate solutions of equations (5) and (6), in conjunction with (7), are given in 
Appendix III. 


3. Initial Conditions 


Strictly, the solutions of equations (5) and (6) must satisfy the condition that 
Q, w, Wand @ vanish when x’,= —°0. However, resort to numerical analysis makes 
it necessary to start the calculations where the induced velocity components U, V (see 
equations (43)) are very small, but not necessarily negligible. Now the analysis of 
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Appendix II shows that, as x’ increases negatively, U tends to zero more rapidly 
than V. Hence, choosing the initial value of A, in Appendix III to correspond with 
some point where U is effectively zero, we have initially 


U,;=0, (4=V;/U, (Vi—> 0), 

and further o,=0, y (9) 
since » <1 in the region concerned. Consequently, from (7), (52) and (53) 

A y=Ay=1, . (10) 


Thus, substituting (8), (9), (10) in (54), setting j=0 in the first interval, and noting 
that A, =0, we obtain 


for0<A<A,. In the case of 9, equation (55) becomes indeterminate, for clearly, 
when A=A,, L,=L,=00. However, if do/dA is regarded as constant in the interval 
concerned—an entirely legitimate assumption—then 


(12) 


Further, combination of (6) and (11) gives 


d 


But, in the region of flow under discussion, (2/2, differs only slightly from unity, and 
sin(y—9)—~(¥—¢). Hence, substituting (13) in (12), we obtain, finally, 


(14) 


in the interval0 <A <A,. 


A value of x’, which satisfies the foregoing argument is x’,= —10c. Equations 
(11) and (14) are then certainly valid for A,=5, ie., for a displacement of the 
aerofoil of Sc. 


4. Impact Velocity and Chordwise Extent of the Roughness 


When 2) and y are known, the development of » and ¢ for individual particles 
may be determined from equations (54) and (55) in combination with (5), (6), (7), 
(44), (45) and (51), (52), (53). Hence, by following this process to a point of impact 
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(if it occurs), and repeating the procedure for a number of particles located initially 
along O’y’, the distribution of , ¢ on the contour of the aerofoil is obtained. 
Clearly, the condition that a particle shall strike the aerofoil may be defined as 


x 
“cos a— ='sin a= Xo J da 
c c c 


x; 
¢ 


y 

sin 008 + | sing da, 
0 


where x, y, are the co-ordinates of a point of impact relative to the aerofoil axes Ox, 
Oy (Fig. 1). Thus, equations (15) determine when and where a particle, initially at 
x’,, Yo, Will make contact. 


Again, relative to the axes Ox,, Oy, of the aerofoil (Fig. 2), the particle velocity 
components are —U,+wcos®# and wsin9, with w<U,. Hence, the resultant 
velocity Wz, in non-dimensional form, is 


Wp =We/U,= (1 — 20 cos ¢ + (16) 


and the tangent to the particle path makes an angle 


_ 


cos ay 


o=tan-' ( 


with Ox,. Thus, the impact velocity w, at x,, vy, has the non-dimensional value 
=O» (i) sin|z + (i) + 7 (18) 


where 2 is the angle of incidence, - the slope of the profile and (i) a cipher to 
indicate the value of a function at the point concerned. The modulus of the angle 
2+ (i)+7 (i) is also taken so that w, remains positive for both surfaces. 


Finally, as already observed in Ref. 3, roughening will cease when w,=(W,)min, 
the latter quantity being the least value of insect rupture velocity relevant to a 


particular flight. (*See footnote on p. 178). Hence, writing “r” in place of “i” to 
indicate impact at the limit of roughness, we have, from (18), 


(«,)min =p sin}z+o (r)+7(r)|. . 
177 
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It remains to relate equation (19) to the streamwise extent of the roughness. 
Due to the difficulty of calculating the induced velocity field set up by an aerofoil of 
arbitrary section, a general approach to the immediate problem is not possible. 
However, the flow about the Joukowski type of aerofoil is easily established 
(Appendix II), and therefore, the corresponding values of x,, y, follow from 
Section 2, in conjunction with Appendix III, and the preceding part of the present 
section, Transition to an aerofoil of any prescribed shape can then be made on the 
assumption that only small differences of y are involved—a condition which is 
sufficiently satisfied with any aerofoil commonly used in practice. Further, from 
(16) and the fact that w < 1, it is evident that, irrespective of aerofoil shape, w, does 
not differ greatly from unity. Consequently, variations of , with respect to both 
x and y may be neglected. Finally, particle paths will remain essentially parallel in 
an area of the flow of order (4x,)’, where 4x,=x,—x,,; and the suffix J refers to the 
basic Joukowski profile. Hence, we may write approximately 


(xs) + (3). 
rd 


7 (x,.) and (dr/dx) ,, being the slope of the arbitrary section and its rate of change 
at x=x,,. Thus, substituting (20) in (19), we obtain, at the limit of roughness for 
any normal shape of aerofoil, 


min = Wry (r) sin|2 + a; (r) (x5) (5) 


or writing, for the moment, 
(xX,3) (r) + 
and noting that, in general, 


sin {a+o(r)+r(r)} ~ 2+o(r)+ r(x), 
(© = Wry (r)\2 + Wy (r)+ Ty + Wpy (r) dr + (=). éx,|. (21) 
rJ 


But the first term on the right hand side of (21) is itself equal to (©,)min. 


to the absolute minimum for a fully representative sample of the aerial population, since clearly 
many random circumstances may arise to cause significant deviations. Nevertheless, in the 
present analysis, it is assumed that, when conditions reach a point at which the roughness 
generated is sufficient to create a serious disturbance to the surface flow, the assumption is 
acceptable, i.e., the population is completely represented and sampled to all intents and 
purposes. Furthermore, it would appear. ® that local differences in the composition of the 
population may be ignored. 
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Hence +f( ix, =0, 


and therefore /dx):,  (dr/dx)z,, 


so that, finally, if x,, is the chordwise extent of the roughness measured positively 
from the leading edge, 


a 7 (ts) 75 (0) 


5. Determination of ©, and ¢; (i) 


The most suitable Joukowski approximation to any general form of aerofoil 
which may be contemplated would seem to be the symmetrical section (6=0), in 
which case the independent variables reduce to a, b, z and «x. Accordingly, the 
tabulation of «,; and oc, (i) from the analyses of Sections 2 and 4 has been under- 
taken in sufficient detail to provide a means whereby approximate solutions for any 
normal association of the relevant parameters can be easily obtained by 
interpolation. 


The relevant calculations were performed for arbitrary values of z equal to 
0°. 5° and 10° and « equal to 0:75 and 3-0, in conjunction with the following 
combinations of a and b: (i) a/c=0-250, b/c=0, (ii) a/c=0-249, b/c=0-015 and 
(iii) a/c =0:247, b/c =0-031, corresponding, respectively, to thickness/chord ratios 
of zero (flat plate), 0:0744 and 0-1540. 


So far as « is concerned, estimates of the quantity k (Fig. 3) from the drag data 
of Ref. 5 indicate that its value near the ground, where roughening takes place 
essentially*’, lies between 0-05 and 0:2 per ft., the lower value referring to relatively 
large insects about 8 mm. long, and the other to the small insects from 1 to 2 mm. 
in length. However, it has been shown'-* that only the small insects are of 
importance in the present problem. Consequently, it would seem reasonable to 
confine & to values in the region of 0:15 to 0-2 per ft. Hence, since the mean wing 
chord commonly lies between 5 and 15 ft., the practical limits of «, i.e. kc, may be 
taken as 0°:75<«<3-0. Further, calculation shows that this parameter is of 
secondary importance, and that values of ,; and oc; (i) for intermediate values of « 
may be interpolated linearly with sufficient accuracy. 


Table I is an abridged record of the data. 


6. Comparison of Theoretical and Experimental Values of x,, 


Careful measurements of x,, on the wings of a Meteor, a Comet and the 
Armstrong-Whitworth E 9/44 (A.W. 52) are recorded in Ref. 2, the observations 
being made after flight. Fig. 4 is a typical example of the way in which x,, develops 
with incidence, and illustrates clearly the feature referred to in Ref. 3, namely, the 
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TABLE I 
VALUES OF 


a=0°. Upper and Lower Surfaces 


Thickness 


| | 
pita 0744 | 
Chord | | 
@-xv/e | «=075 | 300 | O75 | 300 | O75 | 300 
0 0992 | +0953 | 0973 0-921 
0005 0 0 0628 0610 0847 0-802 
0-01 0 0 0:466 0-450 0-733 | 0-698 
0-03 | 0 , 0263 | 0-251 0-486 0-448 
0-05 0 0 0:191 0:176 0-350 0-313 
0-07 0 | 0 0-142 0130 0-224 
0:10 0 0 0-098 0-082 0180 0-143 
ois | | | 0048 | 0-033 0-095 0-058 
020 | oo8 | 0035 | 0-001 
a=5°. Upper Surface 
00 0 0-967 0889, 0946 
0005 0 0-655 0541 0-781 0-697 
0-01 0 0 0368 0269 0-655 
0:03 0 0 0136 | 0071  0:366 0-251 
0-05 0 0 0-034 | 0-219 0-099 
0:07 0 0 | | 0130 0-014 
0:10 0 0 | 0-039 
a=10°. Upper Surface 
0:00 | 0 | 0 | 0955 | 0863 | 0934 0-854 
0-005 0 | 0 0584 | 0477 0-613 
0-01 0 0 0-236 0130 0-447 
0-02 0 0 0-015 0374 0-210 
0-03 0 0 0243 0-057 
0-05 0 0 
a=5°. Lower Surface 
0:00 0130 | 0-198 0-967 0:889 0:946 0866 
0-01 0-122 0-184 0:596 0-602 0-757 0-723 
0-05 0-106 0-151 0-281 0:290 0-442 0-426 
0-10 0-099 | 0-132 0-196 0206 | 0274 0-264 
0-20 0-096 0-114 0-118 0-123 0°136 
0:30 0-094 0-109 0-075 0-078 0-058 0-050 
0-40 0-093 0-106 0050 | 0-054 0-014 0-010 
a=10°. Lower Surface 
0-00 | 0245 | 0336 | 0-955 0-863 0-934 0-854 
0-01 0235 | 0-323 0-655 0-630 0°766 
0-05 0217 | 0-290 0-385 0413 | 0°528 0-530 
0-10 0-206 0:264 0-290 0-315 0-373 0-384 
0-20 0196 | 0-229 0-212 0-228 0-233 0-250 
0:30 0199 | 0209 | 0168 | 0179 0-153 
0-40 | | O199 | ©0142 | | 0-107 0-122 
0°50 O185 | 0182 | O128 | O135 | 0-084 | 0-095 
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TABLE I—continued 
VALUES OF ¢,(i) IN DEGREES 


a=0°. ‘Upper and Lower Surfaces* 


Thickness 0 0-0744 0:1540 
Chord 
(x,-x/e  «=075 | 3:00 075 | 300 | O75 | 3400 
0-00 | o0 | 000 | 000 , 000 
0-005 0 0 016 | 067 035 410 
0-01 0 0 0:20 0-78 052 
0-03 0 0 0:27 1-01 O71 | 2:39 
0-05 0 0 0-31 Vl 0-79 2:70 
0-07 0 0 0-31 1:20 0°84 2:86 
0:10 0 0 0:32 126 0-90 2:98 
0-15 0 0 0:33 1:30 0-92 3-04 
0-20 0 0 0-31 091 | 2286 
a= 5°. Upper Surface 
0-00 0 | 0 2-15 $73 | 57 | $07 
0-005 0 0 2:31 630 , 189 | 718 
0-01 0 0 2:56 6°63 2:30 8°75 
0-03 0 0 3:39 7:17 3:27 10:77 
0-05 0 0 4:10 366 | 10-71 
0-07 0 0 388 10:47 
0-10 0 0 4-06 
a=10°. Upper Surface 
000 380 124 | 896 
0-005 0 0 4:22 11-47 3:33 10-20 
0-01 0 0 5-05 12:64 3°65 12-06 
0-02 0 0 9:17 451 14:89 
0-03 0 0 5-48 16:85 
0-05 0 0 7-42 
a= 5°. Lower Surface 
oo 253 | 723 | 25 | 573 | $07 
0-01 629 164 4-77 0:90 2:49 
0-05 125 | 433 o7 | 227 | 028 0-79 
0:10 085 311 0:50 139 O14 0:33 
0-20 0-64 2:03 0°33 0-98 0-04 ~0:10 
0:30 0-52 1-69 0-23 0°80 -~0-01 
0-40 0-44 1:45 0-19 0-73 ~0-05 ~0-33 
a=10°. Lower Surface 
0-00 525 | 1456 | 3-80 10-24 | 311 | 8% 
0-01 461 | 13-48 3-32 8-41 266 | 719 
0-05 3-47 10-69 2:46 6-42 189 5°55 
0-10 2:75 8°45 1:95 5-23 1:44 4:55 
0:20 2-02 5-70 1-45 3-80 104 | 3-32 
0:30 1:65 4-08 1:15 2:88 083 | 2-70 
0-40 1:36 2-90 1-02 222 | 0-74 2:26 
0°50 | se | | | OF 1:97 


*Note: 7 ,(i) is positive on the upper surface and negative on the lower. 
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Ficures 5 and 6. Comparison of theoretical values of x,, with observations on the Meteor 
(Fig. 5) and A.W.52 (Fig. 6) from Ref. 2. 


rapid increase of x,, on the lower surface, and reverse condition on the upper 
surface, as the incidence increases. Consequently, for the lower surface, the results 
of Ref. 2 indicate the extent of the roughness attained near the minimum speed at 
which the aircraft was flown during the period when roughening took place, and 
conversely, the readings for the upper surface refer to conditions at about the 
maximum speed reached during the same interval. 


Table II gives the relevant experimental data, together with the estimated 
values of x,,, while Figs. 5 and 6 compare theory and observation for the Meteor 
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TABLE II 
SUMMARY OF CALCULATIONS OF X,,. IN RELATION TO THE FLIGHT DATA OF REF. 7 

Aircraft | wing a (x,,)/e 
| surface thickness (m.p.h.) . | 
Meteor Upper | 290 0-083 1:20 0-113 
Lower 0:090c 90 120 0-200 120 0-278 
Upper 0:105¢ 13 290 0-083 146 0-124 
Lower 120 0-200 146 0-287 
Upper 0120e | 13 | 290 0-083 1:75 | 
q Lower 0'120e 90 | 120 | 1:75 | 0-308 
A.W.5S2 Upper | | 250 | 0:096 197 0-211 
Lower O1sse | 77 | 9 0°267 197 | 0-178 
Upper | | 250 | | 2-44 | 0-183 
Lower 0165¢ 87 90 | 0267 | 244 | 0-217 
Upper 0175¢ 08 250 0-096 292 0136 
‘ Lower O17 | 97 | 9 | 0267 | 292 | 0-255 
Comet Upper 2800086 2-73 (0062 
Lower 0-115¢ 74 120 | 0-200 2:73 | 0-192 


and A.W. 52 respectively. Considering the difficulty of obtaining a sharply defined 
measurement of x,,, the agreement in both cases is thought to be as good as can be 
expected. Strictly, the flight conditions refer to three-dimensional flow, particularly 
towards the wing tip, whereas the calculations are confined to two-dimensional flow. 
However, from the experimental and theoretical evidence advanced in Ref. 4, it 
would appear that the effect of a spanwise component of the flow on x,, may be 
ignored unless it is of the order U,, a condition which clearly does not apply to the 
records of Ref. 2. 


The observations made on the Comet were not so complete, and only two values 
of x,. are given, namely, a figure of 0:07c on the upper surface and 0-22c on the 
lower surface, corresponding to a section of the wing at roughly 0-6 of the semi-span. 
Further, the collection of insects was made during the eight stages of a normal 
operational flight trial to Malta, and the speeds relating to the extremes of x,, on 
each surface were not recorded. From subsequent information kindly provided by 
Mr. Johnson, it has been estimated, however, that the result for the upper surface 
relates to a speed between 280 and 300 m.p.h., and, for the lower surface, to a speed 
of about 120 m.p.h. Under these conditions, calculations of x,, yield the results 
given in Table II. Theory and observation, it will be seen, are again in reasonable 
agreement. 


7. Effective Excrescence Height 


The investigations of Ref. 3 showed that the chordwise roughness pattern may 
be sharply divided into two parts, namely, a narrow region adjacent to the leading 
edge where the deposits are highly concentrated and relatively large, followed by a 
considerably more extensive area in which impacts are less frequent and consist 
mainly of fluid traces. It is also evident that the excrescences round the leading 
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edge are invariably of sufficient height to cause premature transition within the 
full-scale range of Reynolds number. The remainder of the roughness proves, 
however, to be insignificant in that, if the relatively restricted section upstream is 
eliminated, or its effects are suppressed, as discussed in Section 5 of Ref. 4, an 
incident laminar layer persists undisturbed. It is of interest, therefore, to establish 
the critical point between significant and insignificant roughness. In this respect, 
the conditions for premature transition due to natural roughness follow essentially 
those associated with artificial roughness of a three-dimensional character and, 
accordingly, depend on the roughness Reynolds number <u./v, where u, is the 
velocity at the top of an excrescence of height «. At the same time, the significant, 
or effective, value of < for the naturally roughened surface is found to be that defined 
by the envelope to the high points. It represents, it would appear, the limiting dis- 
tribution of excrescence height as the result of an indefinitely large number of 
impacts and is satisfied, of course, only at certain points in practice. This feature is 
dealt with more fully in Ref. 4. For the present we are concerned with the condition 
referred to, namely, with the effective value of excrescence height* at which the 
roughness ceases to be a potential source of disturbance to a laminar boundary layer. 


| 


FiGurE 7. Model particle during impact. (a) Instant of contact, w,=w,. (b) Instant when 
w,, =0 for w, < w, (no rupture). (c) Instant when w, =0 for w, > w, (rupture). 


It need hardly be said that impact in the present context is of a very complex 
character, involving the anatomical structure of the insect as well as the purely 
dynamical elements of the phenomenon. At present only a _ very limited 
knowledge exists in regard to the relevant physical factors (see, for example, Ref. 7), 
and a theoretical treatment which takes into account the many diverse elements 
associated with the problem would seem to be incapable of formulation without 
prior research into a number of details. Fortunately, however, the downstream 
roughness, with which we are concerned immediately, appears to be rather simpler 
than that close to the leading edge. As already observed, deposits in the former 
region are composed predominantly of fluid matter, which suggests that the insect 
does not disintegrate, but rather the exoskeletal frame of the body is fractured 
sufficiently to cause a discharge of the haemolymph filling the body cavity. This 
somewhat crude conception of impact has led to the following tentative approach to 
the prediction of effective excrescence height. 


*Since only the effective excrescence height is involved in the subsequent discussion, we shall 
continue to represent it by <. 
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Consider a particle composed of an outer covering (1), an inner membrane (2) 
surrounding a cavity (3), and finally, fluid (4) in the intervening space.* Suppose, 
now, the particle strikes a surface which may be considered inelastic. At the instant 
of contact (Fig. 7(a)), let the mass centre M be a distance n, from the wall, and let 
the velocity of this point in the direction On be w,, with initial value (at contact) 
equal to w;. Then assuming the particle compresses elasticallyt, i.e., the force of 
impact is proportional to the displacement nm of the mass centre, we have 
approximately 


m, being the mass of the insect and c, a constant. Hence, integrating equation (23) 
to satisfy the end conditions n=n,, w,=W, and n=n,, w,=0, we obtain 


where d=n, —n,; and similarly, if d=d, in the particular case when the initial value 
of w, is equal to (W,)min, 


d=(™) 


Thus, by definition, the particle ruptures when d > d,. From previous observations, 
it is now assumed that excrescences are formed from the fluid (4), the volume of a 
deposit, to a first approximation, being proportional to d—d, and the excrescence 
height «. Accordingly, 


and therefore, from equations (18) and (19), 


in which sin {2+ (i, r)+ (i, r)} remains approximately equal to z+ (i, r) + (i, 
but the sign of the angle must now be taken into account. Further, as noted in 
Section 4, w, never differs greatly from unity, so that 
*This model, it will be noted, is a crude representation of the insect components which are 
likely to be of main significance in the present problem; that is to say, (1) corresponds to the 
body shell (here regarded as elastic), (2) and (3) represent the major internal cavity (alimentary 
canal), while the fluid (4) is equivalent to the blood, or haemolymph, which, in circulation, 
fills the entire unoccupied space surrounding the body organs, This amounts to between 15 and 
75 per cent of the total insect volume. 


tEvidence which lends support to this assumption is summarised in Ref. 7, pp. 119-123. 


tThis equation holds strictly for a rigid particle brought to rest by an external spring of 
negligible mass, but has been regarded as an acceptable approximation in the present insiance. 
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and, on condition that x —x,, is small, it is sufficient to write 


()+ (4). (x=), 


r@=r (+ (32). (x-x,). 


Hence, from (24), (25) and (26), we obtain finally 


=constant x (@,)min f (4) (4 — 


do dr 
* 
with f(x)= a+o(r)+7(r) 


For the conditions considered in Ref. 3, it is also evident that the envelope of < 
is an approximately linear function of x in the region with which we are concerned 
and therefore 


Consequently, equation (28) provides an experimental check of the preceding 
analysis. As presented, however, the result expressed by (27) refers to the normal 
case of heterogeneous trapping during flight. Nevertheless, in so far as the argument 
followed is a fair representation of the conditions at impact, it is clear that (27) will 
continue to hold for roughness produced by a particular kind of insect, provided that 
(W,)min is replaced by the appropriate value of w,. Thus, in relation to the results of 
Ref. 3, which were obtained with the fruit fly, Drosophila, it is necessary to take the 
relevant rupture velocity, namely, (w,)p, say. The observations were also made at one 
value of U, (265 ft./sec.), so that (w,)) was constant. Hence, under these conditions, 
f (x,) should be proportional to the measured values of (d:/dx),. As Fig. 8 shows, 
the agreement is surprisingly good, considering the tentative character of the 
assumptions on which equation (27) depends. 


Since the data of Table I refer to free-air conditions. they cannot be used to 
determine f(x,) in the experiments under discussion. It was necessary, therefore, 
to estimate the relevant values of , (i), « (i) from separate step-by-step calculations 
of the kind discussed in Appendix III. These results, together with the measure- 
ments of x,,, also lead to estimates of (w,)) through the equivalent of (19), namely, 


(29) 


(©,)p = sin|a + o (r) + (r)}. 
187 


(26) 
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Further, from Ref. 6, we have a direct measurement of (w,)p but, owing to the 
accelerated state of the insects (see Ref. 3, Section 3), the appropriate value of U, 
cannot be precisely defined. However, a reasonable approximation is 


(w,) 


where U, denotes the tunnel wind speed. Comparisons of (,),, from (29) and (30) 
are given in Table III. 


Considering the difficulty of measuring x,, with precision, and hence the possible 
error in (,)p from (29), due, principally, to the rapid change of 7 over the front part 
of the aerofoil, the agreement is reasonably good, except, perhaps, for N.A.C.A. 
66-009 at small angles of incidence. The discrepancy is thought to be due to 
skidding of the insect on impact, i.e., to a tendency of the insect, when slightly 
ruptured, to slide along the surface, leaving a trace in its path for a short distance. 
Such evidence, from elongation of the deposits, was often noted, particularly when 
the incidence and surface curvature at the downstream limit of roughness were 
small, and similar observations have been made by Johnson”? in flight. It is 
possible, therefore, that only the upstream ends of these particular streaks indicate 


0-07 
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3 
Cf (Xr) 
Ficure 8. Relation of f (x,) to (de/dx), as obtained from the measurements of Ref. 3. 
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TABLE Ill 
___ COMPARISON OF (»,)p) FROM EQUATIONS (29) AND (30) 


Aerofoil Incidence 


From (29) 


Flat Plate 
N.A.C.A. 66-009 


N.A.C.A. 66,018 
| 


points of impact, so that, when they occur at the limit of roughness, their length, by 
the definition of x,,, should be ignored. In fact, however, all measurements of x,, 
were taken to the limit at which the surface became devoid of any sign of contamina- 
tion. Thus, in some instances, x,, may have been somewhat over-estimated, in 
which case 7 in (29), and consequently the corresponding value of w,, will be 
too small. 


The average of (,), from (29), in Table III, is 0-205, or 0-227 at an incidence 
of 6° only, when the most consistent results were obtained. The mean from (30), 
on the other hand, is 0:231. Hence, taking, in general, a figure of 0-225 for (,)p 
and noting, from Fig. 8, that (d¢/dx),/f(x-)=0-0121 in., we obtain, finally, a value 
of 0-054 in. for the constant in (27). 


In Ref. 4, equations (22) and (27) have been applied to show the proportions of 
significant and insignificant roughness under typical flight conditions, the criterion 
for premature transition being based on the relevant analysis of Ref. 3. The 
conclusions reached are in good agreement with those from flight records’, i.e., on 
the upper surface, a maximum of about 5 per cent of the chord from the leading 
edge may be seriously affected by insects, with roughly double this amount on the 
lower surface. 
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Appendix I 


SIGNIFICANCE OF THE FREE ACTIVITY (RANDOM MOTION) OF THE INSECT 


In Section 2 of the main paper, the forces exerted by the insect through wing action 
are neglected, apart from a lift component sufficient to react gravity. Equations (5) and 
(6) then define the hodograph of such a particle. If, however, a propulsive force 7, 
coaxial with the drag D, is generated, the value of « must be based on the resultant force 
D-T. 
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Thus, the power output of the insect is given by 


E—E,=Tq 


where E is the total power and E, the amount required to support the weight of the 
insect. Then it is easily seen that equations (5) and (6) continue to define the hodograph, 
provided that « is replaced by 


kmq 


The zone in which the particle is partially drifting, i.e., the region for which 
(E—E,)/(km,q’) = 1 will now be determined. Measurements of E and E, are given in 
Ref. 5 for a number of different insects and conditions of flight, while the corresponding 
values of k may be determined from Fig. 3 of the present paper. Values of (E—E,)/(km,) 
are found to vary widely but, for the small insects considered here (see Section 5), the 
range is narrow, and a figure of 400 (ft./sec.)® at the maximum speed for prolonged 
insect flight is typical. This condition, it is thought, represents the most frequent 
upper limit of power output which an insect can sustain in all but very brief aerial 
activity near the ground. The calculations may also be simplified by putting g=QU, in 
(31), since » is small compared with 22, and terms containing © can be legitimately 
neglected in the present demonstration. Further, it is sufficient to calculate ( for the 
Joukowski aerofoil (Appendix If), and to ignore the small displacements indicated by 
equations (1). Hence, we take x,/c=—(x’,/c+A) and y,=y’=y’,, the constant y’, being 
chosen to correspond with impact at the limit of roughness, i.e., y’, =x, sin 2+ y, Cos 2. 


o5 


Figure 9. Modification of « due to activity of small insects on (a) the upper surface at 300 m.p.h. 
and 600 m.p.h. (2 ~ 0°) and on (b) the lower surface at 300 m.p.h. (2 ~~ 5°) and 125 m.p.h. 
(z = 10°). 


Figure 9 gives «’/x, as a function of x,, relevant to conditions on the upper and 
lower surfaces within the range of speeds for which roughening is likely to be most 
severe (see, further, Section 2 (a) of Ref. 4). The curves show clearly that, with respect 
to the lower surface (b), small insects begin to act as inert particles when the induced 
velocity is still no more than about two per cent of U, (x,/c=3), and certainly behave 
as such in the important region adjacent to the aerofoil. 


As regards the upper surface (a), it would appear that there is a greater resistance 
to the passage of the aerofoil. However, this is not the case, since, at the small angles 
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Ficure 10. 


of incidence corresponding to the higher speed conditions, induced velocities greater than 
one per cent of U, do not occur at more than about one chord forward of the leading 
edge, i.e., the field is much contracted and therefore the region in which the insect can 
exert a significant influence dynamically is again unimportant. 


Hence, from the foregoing analysis, it is reasonably certain that the power of small 
insects to resist passive drift in the induced velocity field is principally confined to the 
insignificant part of the flow at relatively large distances from the aerofoil. Consequently, 
the assumption of Section 2 that the free activity of insects can be neglected, particularly 
of the small and numerous kinds, is evidently justified. 


Appendix II 
VELOCITY FIELD OF THE JOUKOWSKI AEROFOIL 


For a Joukowski section at incidence z in a stream of velocity U,, parallel to, and in 
the negative direction of, O’ x’, the transformation relations between the aerofoil in the 
z-plane and the circle of radius R in the (-plane (Figs. | and 10) are 


2 

C(=f£+in=re® e ° e e (34) 


in which, from Fig. |, the co-ordinates of P are 


(35) 


x= (x’—x',—U,1) cos 2+’ sin 


y= sin cos , 
x’, y’ being given by equations (1). Then, combination of equations (33) and (34) yields 
464 — + (Sx? + y? + 4a’) (x? . (36) 


and n= = ‘ . en 
as equations for £, » in terms of x, y. 
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There are also the co-ordinate relations (Fig. 10) 


where f,=& cos 2—n sin z—b cos (2 +8) 
(40) 


f,=£ sin 2 +n cos 2—bsin(z+8) . 


Hence, the induced velocity components u (z,), v (z,) in the directions Ox,, Oy, are given 
by the velocity function 


dF dF di, d& dz 


u (z,)—iv (z,)= dz, = de dz, ’ (41) 
dF iK 
with 


and the circulation K, by Joukowski’s hypothesis, equal to 4xRU,, sin (z+ 3), so that, 
finally, we obtain from (41), in conjunction with (33), (38) and (39), 


q,°+4q,? 

__ Pi P2 
U, +4," 


(42) 


where Pp, h, Enh, 
Pa = hy 
{f.+R sin(z+8)} 
h,=2f, {f.+R sin(z+)} . 

Thus, the induced velocity components U, V are 


(43) 


Appendix III 
APPROXIMATE SOLUTIONS OF EQUATIONS (5) AND (6) 


To obtain the distributions of w and ¢ on the contour of the aerofoil, and hence the 
impact conditions, it is necessary to solve equations (5) and (6) step by step from known 
initial values (see Section 3). This may be performed conveniently by iteration which, 
for a system with n degrees of freedom, is defined by 


dQ, 
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in the interval (A—A,), Q,'' being arbitrary. However, where a considerable number of 
specific conditions require to be investigated, as in the present instance, it is particularly 
desirable that the first approximation shall, if possible, be sufficiently accurate to make 
the iteration process a relatively small, or negligible, part of the computation. From this 
point of view, the form in which the hodograph equations have been expressed proves to 
be the most suitable, for, outside a restricted region in the vicinity of the leading edge, 
where /, changes rapidly, the following linear relations are found to hold with good 
accuracy, namely, 


(47) 


a,, b, being treated as constants in the interval (A—A,), and the prime denoting 
differentiation with respect to A. Hence, substituting (46), (47), (48) in (5) and (6), we 
obtain. to the accuracy of first order differences in A,, A,, 


a? {(l—a,) — do=K{A,, +A, dA. . (49) 


i—b, sin(¥—o)~ LA,, «A? (Q,—0,)? 


Further. we have generally 
while we obtain from (7) 
LO (Q+o0) QA, A, @A, 
A, 


Te (Q+) OA, _ A, 
~ A, (Q—e)*’ en 


__ OA, _ A, oA, 
A, (Q—w)? a 6A, OF 


and from (7), (46) and (47) 


Equations (49) and (50) may be solved by quadrature for the values of a,, 6, which, 
when substituted in (46) and (47), finally yield 


KO (Q)—w,) (A,,+4A,, (A—A,)] (—A) (54) 
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L, 
55 
Ay Ags 0) 
AN—Ao 
A, A, ANAw 
yA —A,,A,/ AN=0-9Q,, Aw=o — . 


For moderate values of (A—A,),* equations (54) and (55), in conjunction with (52) 
and (53), lead to values of w and @ which are of sufficient accuracy everywhere other 
than in the region referred to earlier. Near the leading edge, therefore, the foregoing 
solution must normally be regarded as a first approximation to be improved by the 
iteration procedure defined by equations (44) and (45). 


It will also be noted that (55) cannot be solved directly, since L, involves @, but 
cyclical application of (47), (57) and (55) provides successive approximations which are 
rapidly convergent. 
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Some Simple Transformation Functions for 
Square and Triangular Holes with 
Rounded Corners 


W. H. WITTRICK, M.A., Ph.D., F.A.A., A.F.R.Ae.S., A.F.LAS. 


(Department of Aeronautical Engineering, University of Sydney) 


Summary: Transformation functions, in the form of three-term polynomials, 

are derived for the conformal mapping of the region outside a square or an 

equilateral triangle, having rounded corners with varying and controllable 
degrees of curvature, on to the region outside the unit circle. 


1. Introduction 


A problem of considerable interest in aircraft design is the determination of 
the stress concentration arising near holes, such as windows in pressure cabins. 
The problem can be solved by complex variable methods if the region outside the 
hole can be mapped by means of a polynomial transformation function on to the 
region outside the unit circle. 


Shapes of practical interest are the ellipse, for which a simple and exact trans- 
formation is available, and the square and equilateral triangle with rounded corners. 


It is well known that the Schwarz-Christoffel transformation enables the region 
outside any polygon to be mapped on the region outside the unit circle. This 
transformation function is in the form of an indefinite integral, which can be 
expanded into an infinite series. If only the first few terms are used the effect is 
to produce an approximation to the polygon with the corners rounded; the greater 
the number of terms, the smaller does the radius of curvature at the corners 
become. It is apparent, however, that in practice this transformation suffers from 
two great disadvantages : — 


(a) The smallest possible number of terms required to give a reasonable ap- 
proximation to a regular polygon is two, to which corresponds a definite 
corner curvature. If more terms are taken the corner curvature increases. 
Therefore, if a smaller corner curvature is required, the transformation is 
of no use. 


Corresponding to each number of terms taken, there is a definite value 
for the corner curvature and there is no apparent way of producing inter- 
mediate values. 
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Figure 1. Square 
holes with rounded 
corners, showing one 
corner only. 


This paper is essentially concerned with deriving a modified transformation 
function which overcomes these objections in a very simple way. 


2. The Square Hole with Rounded Corners 


Consider an exact square hole in the z-plane, with centre at the origin and sides 
parallel to the real and imaginary axes. The Schwarz-Christoffel transformation, 
which maps the region outside the square on to the region outside the unit circle 
in the (-plane, is 


where the real constant R determines the size of the square. 


The centres of the sides of the square correspond to the points (= +1 and +i 
on the circle, while the corners correspond to (=exp(+iz/4) and exp( + 3iz/4). 


Instead of equation (1), take a transformation function containing three terms 
only, similar to the first three of equation (1) but with, at present, arbitrary 
coefficients. 
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TRANSFORMATION FUNCTIONS 


It is easy to see that, if R, p and q are real, the inclusion of only these powers 
of ¢ gives a shape symmetrical both about the real and imaginary axes and about 
the two lines at 45° to these axes, as required. 


As before, the constant R determines the size of the hole. The constant g 
will be determined as a function of p by prescribing zero curvature at the centre of 
each side of the “square”. The remaining constant p can then be varied to control 
the corner curvature. 


Let @ be the polar angle in the (-plane, so that, on the unit circle, 
(=exp (i6). Substituting into equation (2) and separating real and imaginary parts, 
the parametric equation of the hole in the z-plane becomes 


x= R (cos 6— pcos 36 + q cos 76) 
y= R (sin 6+ p sin 34 — q sin 76) 
where z= x+ ly. 


The curvature 1/r of the hole boundary (positive when concave inwards) is 
readily shown to be 


where a dot denotes differentiation with respect to 4. 


Using equations (3) and (4), the curvature at the points 6=0, = and +7/2, 
which will correspond to the centres of the sides of the “square”, can be shown 
to be 


1 (1-9p+49q) 


R (1+3p—7q)° 
The constant g is now chosen to make this zero, so that 

Again using equations (3) and (4), the curvature 1/r. at the “corners” of the 


“square”, which correspond to the points 6= +7/4 and +3/4, is given by the 
equation 


On using equation (5), this becomes 


r.  2(4-15p)?° 
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Instead of using the ~onstant R to control the size of the hole it is convenient 
to work in terms of the “width” 26 of the “square”, defined as the distance between 
the centres of opposite sides. Alternatively. b may be defined as the radius of the 
inscribed circle. Thus the points (=1 and z=b correspond, so that, from 
equation (2), 


b=R(l1-p+q). 
Using equation (5) this gives 


b 8 


Hence, on using equations (2), (5), (6) and (7), the final transformation function 
becomes 


z 4% 9p-1)i-" 


8 (6—5p) ©) 
while the corner radius of curvature is given by 

_ 


Some members of the family of holes, corresponding to different values of p 
in equation (8), are shown plotted in Fig. 1. It will be seen that the central portions 
of the sides are remarkably straight, and in this respect it is relevant to note that 
not only is the curvature forced to be zero at the centre of each side by virtue of 
equation (5), but also the symmetry automatically ensures that the rate of change 
of curvature there is also zero. 


Note also that one member of the family, for which p = 1/9, corresponds to only 
two terms in the transformation function, its corner radius of curvature being 6/4. 


3. The Equilateral Triangular Hole with Rounded Corners 


An exactly similar procedure can be used to obtain a family of holes in the 
shape of an equilateral triangle with rounded corners. The centre of the triangle 
is at the origin in the z-plane and one side is parallel to the imaginary axis on the 
negative side. The size of the hole is defined by the radius b of the inscribed 
circle of the triangle. 


The final transformation function, which gives zero curvature at the centre of 
each side, is 


Z_ 250+ + 
b 3(8-7p) 
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Ficure 2. Equilateral triangular holes with rounded corners, showing one corner only. 


The radius of curvature r, at the “corners” is given by the equation 


(3 —7p) 


~ 6p 8-0) 


A few members of the family, corresponding to different values of p. are shown 
plotted in Fig. 2. 


Note that again there is one member of the family, with p=1/4, for which the 
transformation function contains only two terms; the corresponding corner radius 
of curvature is 
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REVIEW 


Review 


The Scientific Papers of Sir Geoffrey Taylor. Vol. 1, Mechanics of Solids; Vol. H, 
Meteorology, Oceanography and Turbulent Flow. Edited by G. K. Batchelor. Cambridge 
University Press. 600 pp. and 500 pp. Illustrated. 75s. each. 


Before commencing a new line of investigation a research worker must first find out 
what has already been written on the subject, both by consulting books and, above all, 
by reading the original research papers. An enquiry starts by searching the literature 
and, more often, by asking others for the names of scientists who are likely to have made 
contributions in a particular field. In fact, a lead towards the names is often the best 
way to start. The same is true to a lesser extent of anyone who wants information on a 
given subject. Much of the key work will be found in this way and the bibliographies 
included in the research papers will direct further reading. Thus, in the early days of 
this century, the mathematical physicist often consulted the collected papers of Lord 
Rayleigh. Equally today a physicist would turn to the collected papers of Sir Geoffrey 
Taylor, whose clear exposition and breadth of outlook, in a similar but broader realm, 
have rarely been equalled. 


G. I. Taylor’s collected papers will fill four volumes and these first two, containing 
86 papers together, have been admirably produced as regards text, mathematics, line 
diagrams and photographic reproductions. They have been carefully edited. Our regret 
is the absence of any clue as to what is in each volume other than the undertaking by the 
editor, Dr. G. K. Batchelor, that a complete list of the papers in all four volumes will be 
provided at the end of Volume IV. Sir Geoffrey is the 18th in a most distinguished line 
of scientists whose collected papers have been published by this press. The list of 
Contents in each volume gives the title of the paper, its co-author if any, its date and 
place of original publication. A high proportion were published by the Royal Society, 
a few were published abroad and some are now made generally available for the 
first time. 


Volume I commences with papers on the use of soap films for solving torsion and 
flexure problem and were written while he was stationed at the R.A.E. in the First World 
War. Many papers follow on the distortion and plastic extension of single crystals of 
aluminium, iron, beta brass, etc. The general problems of the shear and distortion 
of metals retained his interest over many years, and after a break during the Second 
World War he returned again to the subjects of elasticity and plasticity. His interest in 
the war effort was shown by papers on the propagation of earth waves in an explosion, 
the mechanical properties of cordite, and the effects of impact loads. His James Forrest 
lecture was entitled “The testing of materials at high rates of loading’. This volume 
contains his famous work on dislocation theory. 


The phenomenon of turbulence and its effects in the atmosphere and in problems of 
fiuid mechanics provides the main theme of the second volume, but it also contains a wide 
range of geophysical problems. The early papers on meteorology indicate the direction 
in which his ideas evolved up to his well-known work on the statistical theory of 
turbulence. The earliest on eddy motion in the atmosphere was published in 1915 and 
was soon followed by others on turbulence in the lower atmosphere and on the dissipa- 
tion of eddies. His five papers on Statistical Theory and that on the Spectrum of 
Turbulence, published by the Royal Society between 1935 and 1938, are all included in 
this volume. We find also papers on diffusion in the atmosphere, formation of fog, tidal 
friction and tidal oscillations, and the stability of stratified fluids. In fact, this volume 
will interest meteorologists, oceanographers, and geophysicists, as well as scientists 
concerned with fluid mechanics. J. L. NAYLER. 
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Aircraft 


EXPANDING COMMITMENTS LEAD TO AN: 
_ EXPANDING GUIDED WEAPON TEAM 


ADVANCED WORK ON 
HYPERSONIC PROJECTS 


AERODYNAMICS— 


Aeronautical Engineers or 
Mathematicians 


Al. Advanced problems of hypersonic 
flow and their application to configuration 
design. 


A2. Unsteady flow aerodynamics. 

A3. Stability and response. 

A4. Supersonic project aerodynamics. 
A5. Hypersonic flutter. 


STRUCTURES— 
Engineers or Math 


Sl. Project design of advanced missiles. 


$2. Heat transfer problems. 


INSTRU MENTS— 
Physicists 
Pl. Optical devices, including mosaic 
systems, colour densitometry and flight 
camera installations. 
P2. Gyroscopes, pressure measuring de- 
vices and electro-mechanical instruments. 
FLIGHT SYSTEMS— 
Fl. Mathematicians or Engineers for 


theoretical studies and development on 
missile servo-control systems. 


F2. Mechanical Engineers for theoret- 
ical studies and development on hydraulic 


and fuel systems. 
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F3. Physicists or Electronic Engineers 
for radome development. 


F4. Physicists or Electronic Engineers 
for aerial development. 

ELECTRONIC TECHNIQUES— 
Physicists 
El. Development of techniques associated 
with advanced electronic devices. 


WEAPON SYSTEMS— 

Mechanical and Electrical Engineers 

WI. Liaison with manufacturers of the 

Bloodhound Weapon System in connection 

with overseas sales. 
These extremely interesting opportunities are 
offered by a member of Europe’s largest 
guided weapon consortium, which has 
recently received large orders for Bloodhound. 

Applicants should have a degree or similar 
qualifications and 2 or 3 years’ experience in 
related fields. Salaries will match the high 
standards required and assistance will be 
given with housing and removal expenses. 

Please write to Mr J. Raimes, giving full 
details of age, qualifications and experience and 
quote the reference numbers of the posts in which 
you are interested. Further information on these 
posts will be supplied on request. 

J. Raimes, Personne! Manager, 

Bristol Aircraft Ltd, Room AEO/60/AQ, 

Filton, Bristol. 
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if you are an A.F.R.Ae.S. or possess an equivalent 
qualification, please write to the Technical Staff Manager, 
when we shall be pleased to give you details of 
the opportunities available to experienced men within our 
technical organisation. 


Blackburn Aircraft Limited, Brough, Yorkshire 
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50 years of rience. 
y expe 
f 
are behind the Blackburn NA,39 
1908 Robert Blackburn designed, built and flew his first. 
aircraft at Marske Sands Yorkshire. 
bat 
1959 Latest of a long line, the Blackburn NA.39 Strike Ai ratt,. 
Of this aireraft the Miniatey of Defence has 
“tm the low-level strike role, the of any 


V-bombers fly over 10 miles a minute, 


over 10 miles up, 


for over 10 hours... 


...AND BRISTOL SIDDELEY 
SUPPLY THE POWER 


One of the largest manufacturers of motive 
power units in the world, Bristol Siddeley 
Engines Limited produce two outstanding 
high-thrust turbojet engines—the Olympus 
and the Sapphire. Between them these two 
engines power the major part of the RAF’s 
strategic V-bomber force. The Olympus, 
which delivers 17,000 Ib thrust dry, powers 
the Avro Vulcan Mk 1 and 2. The latest 
version, the Olympus 21, is rated at 20,000 
Ib thrust without reheat. The Sapphire 
powers the Handley Page Victor Mk 1. 
These engines give the V-bombers sonic 
capability — long range — great altitude —- 
superior performance to any other aircraft 
of their type in the world. 
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Bristol Siddeley products currently in produc- 
tion include :— 


AERO 
TURBOJETS: Olympus, Sapphire, Orpheus, 
Viper. 
TURBOPROPs: Proteus, Double Mamba. 
ROCKETS: Gamina. 
RAMJETS: Thor. 


INDUSTRIAL & MARINE 
DIESELS: Maybach. 
GAS TURBINES: Proteus. 
cars: Armstrong Siddeley Star Sapphire, 
Bristol 406, 
OTHER PRODUCTS: Gas bearing compressors, 
Beaver ball screws and splines. 
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ELLIOTT BROTHERS (LONDON) LTD. 


RADAR AND COMMUNICATIONS 


| Automatic Landing 
Systems 


A SENIOR PHYSICIST OR ENGINEER 


is to be appointed in connection with a new project on automatic landing. 

Applicants should have a degree or equivalent together with appropriate 
l experience in the microwave and/or electronics field, and in addition 
should possess the drive, initiative and enthusiasm required to fully 
exploit the project in question. 


Please address replies to:-— 
Personnel Manager (Ref.: 121), 
Elliott Brothers (London) Limited, 
Elstree Way, Borehamwood, Herts. 
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NOTICE TO CONTRIBUTORS 


All communications should be addressed to the Editor, The Aeronautical 
Quarterly, 4 Hamilton Place, London, W.1. 


Papers for consideration should be typewritten with double and be 
accompanied by a summary of not more than 250-300 words. Papers must be 
as short as possible, otherwise considerable delay in publication will result, and 
in any case they should not be longer than twenty foolscap pages. Full 
descriptions and —— results of work, as written for an internal report, are 
not suitable for publication in The Aeronautical Quarterly, and such reports 
need to be rewritten in concise form. Only necessary equations should be 
included in the text; routine mathematics should be omitted or given as an 
Appendix to the paper. Illustrations should be reduced to a minimum. The 
titles of papers should be kept short—with sub-titles if necessary. 


Only simple symbols and formulae should be typewritten. All others should 
be written clearly in ink. Greek letters should be indicated and a clear distinction 
drawn between capital and small letters. Particular attention should be 

to the use of alpha and a, kappa and k, mu and u, nu and »v, eta and n. 
— the numeral 1, and the letter O from the 
zero 0. 


The use of dots, bars, and so on, over letters, or the use of dots as multi- 
plication signs and bars for brackets should be avoided as far as possible. 


Suffixes and prefixes must be indicated clearly and complicated suffixes 
and prefixes should be avoided. 


. A list of all symbols used, both in the illustrations and text, must be given at the 


end of the paper whether they are standard or not. 


. Clear diagrams or sketches must be submitted for redrawing by the Society’s 


draughtsman. 
Photographs must be not less than half plate size and must be clear black 
and white glossy prints. 


. References should be given in the form: — 


Warp, G. N. Calculation of Downwash Behind a Supersonic Wing. The 
Aeronautical Quarterly, Vol. 1, pp. 35-38, May 1949. 


. The copyright of every paper printed in The Aeronautical Quarterly is the 


property of The Royal Aeronautical Society. Permission to reprint or to use 
any paper will not be unreasonably refused. 


. It is a condition of publication that the author shall have obtained consent where 


necessary to use any material in his paper which is copyright or the property 
of any other person or his employers. 


. Authors will be entitled to 20 reprints of their papers free of charge. This 


number will be shared between joint authors. 


at the Offices of The Royal 


Telephone: Grosvenor 3515 (5 lines). Telegraphic Address: Didaskalos, Audley, London 


| 
| 
6 
; 4 Hamilton Place, Piccadilly, W.1, at 15s. Od. net for a single number to non-members 
; and 7s. 6d. net for a single number to members of The Royal Aeronautical Society. Instructions 2 
: for Contributors are given above. None of the papers must be taken as expressing the 
i opinion of the Council of The Royal Aeronautical Society unless such is definitely stated to 
be the case. All communications for publication or advertisement in The Quarterly should be 
addressed to—The Editor, The Acronautical Quarterly, 4 Hamilton Place, W.1. 


ELLIOTT BROTHERS (LONDON) LTD. 


RADAR AND COMMUNICATIONS 


Automatic Landing 


Systems 


A SENIOR PHYSICIST OR ENGINEER 


is to be appointed in connection with a new project on automatic landing. 

Applicants should have a degree or equivalent together with appropriate 

experience in the microwave and/or electronics field, and in addition 

should possess the drive, initiative and enthusiasm required to fully 
exploit the project in question, 


Please address replies to: 
Personnel Manager (Ref.: 121). 
Elliott Brothers (London) Limited, 
Elstree Way, Borehamwood, Herts. 


The Aeronautical Quarterly (Iv) May 1960 


NOTICE TO CONTRIBUTORS 


All communications should be addressed to the Editor, The Aeronautical 
Quarterly, 4 Hamilton Place, London, W.1. 


Papers for consideration should be typewritten with double spacing and be 
accompanied by a summary of not more than 250-300 words. Papers must be 
as short as possible, otherwise considerable delay in publication will result, and 
in any case they should not be longer than twenty foolscap pages. Full 
descriptions and complete results of work, as written for an internal report, are 
not suitable for publication in The Aeronautical Quarterly, and such reports 
need to be rewritten in concise form. Only necessary equations should be 
included in the text; routine mathematics should be omitted or given as an 
Appendix to the paper. Illustrations should be reduced to a minimum. The 
titles of papers should be kept short—with sub-titles if necessary. 


Only simple symbols and formulae should be typewritten. All others should 
be written clearly in ink. Greek letters should be indicated and a clear distinction 
drawn between capital and small letters. Particular attention should be paid 
to the use of alpha and a, kappa and k, mu and u, nu and v, eta and n. The 
poy must be distinguished from the numeral 1, and the letter O from the 
zero 0. 

The use of dots, bars, and so on, over letters, or th: use of dots as multi- 
plication signs and bars for brackets should be avoided as far as possible. 

Suffixes and prefixes must be indicated clearly and complicated suffixes 
and prefixes should be avoided. 


A list of all symbols used, both in the illustrations and text, must be given at the 
end of the paper whether they are standard or not. 


Clear diagrams or sketches must be submitted for redrawing by the Society’s 
draughtsman. 

Photographs must be not less than half plate size and must be clear black 
and white glossy prints. 


References should be given in the form: — 
Warp, G. N. Calculation of Downwash Behind a Supersonic Wing. The 
Aeronautical Quarterly, Vol. 1, pp. 35-38, May 1949. 


The copyright of every paper printed in The Aeronautical Quarterly is the 
property of The Royal Aeronautical Society. Permission to reprint or to use 
any paper will not be unreasonably refused. 


It is a condition of publication that the author shall have obtained consent where 
necessary to use any material in his paper which is copyright or the property 
of any other person or his employers. 


Authors will be entitled to 20 reprints of their papers free of charge. This 
number will be shared between joint authors. 


The Aeronautical Quarterly is published at the Offices of The Royai Aeronautical Soctety, 
4 Hamilton Place, Piccadilly, London, W.1, at 15s. 0d. net for a single number to non-members 
and 7s. 6d. ne¢ for a single number to members of The Royal Aeronautical Society. Instructions 
for Contributors are given above. None of the papers must be taken as expressing the 
opinion of the Council of The Royal Aeronautical Society unless such is definitely stated to 
be the case. All communications for publication or advertisement in The Quarterly should be 
addressed to—The Editor, The Aeronautical Quarterly, 4 Hamilton Place, W.1. 


Telephone: Grosvenor 3515 (5 lines). Telegraphic Address: Didaskalos, Audley, London 


7. 


THE AERONAUTICAL QUARTERLY 


It is the aim of The Aeronautical Quarterly to attract not only original papers contributing 
to aeronautical science and engineering, and papers developing new or improved methods of 
analysis and experimental techniques, but also papers on allied sciences which have a bearing 
on aeronautical problems. The Aeronautical Quarterly is open to authors of any nationality 
and is not restricted to members of the Society. 


Most papers can be published within nine months of receipt. Short papers. of not more 
than 2,000 words, are published much more quickly. 


VOLUME Xl 


Volume XI of The Aeronautical Quarterly will be published in four parts, in February, 
May, August and November 1960. The prices, including postage and packing, are :— 


To Members of the Society To Non-Members 

0.6. 
Per Part 8 0 ($1.15) 15 6 ($2.20) 
Subscription (4 parts) . 112 0 ($4.50) 3 2 ($8.70) 


Papers to be published in forthcoming issues include : — 


C. Bourque and B. G. Newman—Reattachment of a Two-Dimensional, Incompressible 
Jet to an Adjacent Flat Plate. 


W. H. Wittrick—Analysis of Stress Concentrations at Reinforced Holes in Infinite Sheets. 


A selection of the more specialised papers read before the Society during the current 
session will also be published in the Aeronautical Quarterly. 


PRINTED BY THE LEWES PRESS, WIGHTMAN & CO. LTD., LEWES, SUSSEX, ENGLAND, AND PUBLISHED 
BY THE ROYAL ABRONAUTICAL SOCIETY, 4 HAMILTON PLACE. LONDON, W.1, ENGLAND. 


| : 


